For R beginners

New chunk Ctrl+Alt+I

Execute chunk Ctrl+Shift+Enter

Execute all chunks Ctrl+Alt+R

HTML preview Ctrl+Shift+K

Library preparations

library(readr)
library(dplyr)

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union
library(tidyverse)
── Attaching core tidyverse packages ────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse 2.0.0 ──
✔ forcats   1.0.0     ✔ stringr   1.5.1
✔ ggplot2   3.5.1     ✔ tibble    3.2.1
✔ lubridate 1.9.3     ✔ tidyr     1.3.1
✔ purrr     1.0.2     ── Conflicts ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errors
library(ggplot2)
library(reshape2)

Attaching package: ‘reshape2’

The following object is masked from ‘package:tidyr’:

    smiths
library(stats)
library(tm)
Loading required package: NLP

Attaching package: ‘NLP’

The following object is masked from ‘package:ggplot2’:

    annotate
library(text2vec)
Registered S3 method overwritten by 'data.table':
  method           from
  print.data.table     
library(textstem)
Loading required package: koRpus.lang.en
Loading required package: koRpus
Loading required package: sylly
For information on available language packages for 'koRpus', run

  available.koRpus.lang()

and see ?install.koRpus.lang()


Attaching package: ‘koRpus’

The following object is masked from ‘package:tm’:

    readTagged

The following object is masked from ‘package:readr’:

    tokenize
library(syuzhet)

Data Import

data_posts <- read.csv("~/4year/2semester/dtII/CSVs/HEIs.csv",
                 colClasses = c(tweet_id = "character"))

# Modifying created_at type so that attribute can be used more easily 
data_posts$created_at <- as.POSIXct(data_posts$created_at,
                              format= "%Y-%m-%dT%H:%M:%S", tz="UTC")

#View(data)
summary(data_posts)
      id              tweet_id             text               type           bookmark_count    favorite_count     retweet_count      reply_count      
 Length:11728       Length:11728       Length:11728       Length:11728       Min.   :  0.000   Min.   :    0.00   Min.   :   0.00   Min.   :   0.000  
 Class :character   Class :character   Class :character   Class :character   1st Qu.:  0.000   1st Qu.:    7.00   1st Qu.:   2.00   1st Qu.:   0.000  
 Mode  :character   Mode  :character   Mode  :character   Mode  :character   Median :  0.000   Median :   20.00   Median :   5.00   Median :   1.000  
                                                                             Mean   :  1.543   Mean   :   60.67   Mean   :  10.62   Mean   :   3.888  
                                                                             3rd Qu.:  1.000   3rd Qu.:   57.00   3rd Qu.:  11.00   3rd Qu.:   3.000  
                                                                             Max.   :418.000   Max.   :41655.00   Max.   :4214.00   Max.   :2317.000  
                                                                                                                                                      
   view_count        created_at                       hashtags             urls            media_type         media_urls       
 Min.   :      5   Min.   :2022-08-01 03:05:11.00   Length:11728       Length:11728       Length:11728       Length:11728      
 1st Qu.:   2643   1st Qu.:2022-10-19 12:56:27.00   Class :character   Class :character   Class :character   Class :character  
 Median :   6240   Median :2023-01-29 08:26:30.00   Mode  :character   Mode  :character   Mode  :character   Mode  :character  
 Mean   :  14182   Mean   :2023-01-30 07:39:34.96                                                                              
 3rd Qu.:  16058   3rd Qu.:2023-05-05 14:16:43.25                                                                              
 Max.   :7604544   Max.   :2023-08-31 20:50:01.00                                                                              
 NA's   :4840                                                                                                                  

Initial Data Preparation

# Count of how many entries each HEI has
number_posts <- data_posts %>%
              group_by(id) %>% summarise(count = n())

number_posts

Since complutense only has 1 entry we can’t learn anything from it, so we removed it

data_posts <- data_posts[data_posts$id != "complutense.csv", ]

Visualization of all posts, just tweets and just replies

number_posts <- data_posts %>%
              group_by(id) %>% summarise(posts = n())

number_tweets <- data_posts[data_posts$type == "Tweet", ] %>%
              group_by(id) %>% summarise(tweets = n())

number_replies <- data_posts[data_posts$type == "Reply", ] %>%
              group_by(id) %>% summarise(replies = n())

print(number_posts)
print(number_tweets)
print(number_replies)

Calculating the percentage of tweets and replies based on all posts

# Merging the counts of tweets and replies with the count of posts
data_ratio <- merge(number_posts, number_tweets, by = "id", all = TRUE)
data_ratio <- merge(data_ratio, number_replies, by = "id", all = TRUE)


data_ratio$percentage_tweets <- round(((data_ratio$tweets / data_ratio$posts) * 100), 2)
data_ratio$percentage_replies <- round(((data_ratio$replies / data_ratio$posts) * 100), 2)

data_ratio <- data_ratio[, c("id", "percentage_tweets", "percentage_replies")]

data_ratio$percentage_replies[is.na(data_ratio$percentage_replies)] <- 0

print(data_ratio)

NA removal

Function to visualize the number of NAs in all columns

na_count <- function(){
  # Counting the number of NA values for each column
  na_count <- colSums(is.na(data_posts))
  
  # Creating a new data frame with the NA counts
  na_counts_table <- data.frame(Column = names(na_count), NA_Count = na_count)
  
  print(na_counts_table)
}

Calculations of view, favourite, retweet and reply percentiles and visualization of NAs in all columns

data_posts <- data_posts %>%
  group_by(id) %>%
  mutate(view_percentile = ntile(view_count, 100),
         favorite_percentile = ntile(favorite_count, 100),
         retweet_percentile = ntile(retweet_count, 100),
         reply_percentile = ntile(reply_count, 100)) %>%
  rowwise() %>%
  mutate(avg_percentile = round(mean(c(view_percentile, favorite_percentile, retweet_percentile, reply_percentile), na.rm = TRUE), 2))

na_count()

data_percentile <- data_posts[, c("id", "view_percentile", "favorite_percentile", "retweet_percentile", "reply_percentile", "avg_percentile")]

print(data_percentile)

Calculation of the maximum number of views for each HEI

max_view_counts <- tapply(data_posts$view_count, data_posts$id, max, na.rm = TRUE)

print(max_view_counts)
      duke.csv       epfl.csv        goe.csv    harvard.csv  leicester.csv manchester.csv        mit.csv         sb.csv   stanford.csv    trinity.csv 
        307969         105095          15455        2982704          47838         317086        7604544         607498         222593         205333 
        wv.csv       yale.csv 
        109265         143108 

Removal of NAs

# From view count
data_posts$view_count <- ifelse(
  is.na(data_posts$view_count),
  round(max_view_counts[data_posts$id] * (data_posts$avg_percentile / 100)),
  data_posts$view_count)

# From view percentile
data_posts$view_percentile <- ifelse(
  is.na(data_posts$view_percentile),
  data_posts$avg_percentile,
  data_posts$view_percentile)

Visualization of NAs in all columns

na_count()

Function to calculate average posts

average_posts <- function(timeframe){
  # Calculation of the timeframe between earliest and latest post for each HEI
  date_range <- data_posts %>%
    group_by(id) %>%
    summarise(min_date = min(created_at),
              max_date = max(created_at)) %>%
    mutate(num_days = as.numeric(difftime(max_date, min_date, units = timeframe)))
  
  # Naming the column respecting the timeframe
  column_name <- paste0("avg_posts_per_", timeframe)
  
  # Calculation of the number of posts per day for each HEI
  posts_per_timeframe <- number_posts %>%
    left_join(date_range, by = "id") %>%
    mutate(!!column_name := round((posts / num_days), 2))
  
  print(posts_per_timeframe)
  return(posts_per_timeframe)
}
posts_per_day <- average_posts("days")
posts_per_week <- average_posts("weeks")

Plot for the average number of posts per day for each HEI

barplot(posts_per_day$avg_posts_per_days,
        names.arg = posts_per_day$id,
        main = "Average Posts per Day",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(posts_per_day$avg_posts_per_days) + 1),
        las = 2,
        col = "#3498DB")

# Adding text labels over each bar and aligning it with the center of each bar 
text(x = barplot(posts_per_day$avg_posts_per_days, plot = FALSE),
     y = posts_per_day$avg_posts_per_days,
     labels = round(posts_per_day$avg_posts_per_days, 2),
     pos = 3)

Plot for the average number of posts per week for each HEI

barplot(posts_per_week$avg_posts_per_weeks,
        names.arg = posts_per_week$id,
        main = "Average Posts per Week",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(posts_per_week$avg_posts_per_weeks) + 5),
        las = 2,
        col = "#E74C3C")

text(x = barplot(posts_per_week$avg_posts_per_weeks, plot = FALSE),
     y = posts_per_week$avg_posts_per_weeks,
     labels = round(posts_per_week$avg_posts_per_weeks, 2),
     pos = 3)

Defining the intervals of time for the academic year

intervals <- list(
  interval1 = as.POSIXct(c("2022-08-31", "2022-12-15")),
  interval2 = as.POSIXct(c("2023-01-04", "2023-04-01")),
  interval3 = as.POSIXct(c("2023-04-14", "2023-06-15"))
)

Function to check if a date falls within a given interval of time and apply appropriate Boolean

check_interval <- function(date) {
  for (i in 1:length(intervals)) {
    interval_start <- intervals[[i]][1]
    interval_end <- intervals[[i]][2]
    if (date >= interval_start & date <= interval_end) {
      return(TRUE)
    }
  }
  return(FALSE)
}
data_posts$academic_year <- sapply(data_posts$created_at, check_interval)
print(data.frame(id = data_posts$id, academic_year = data_posts$academic_year))

Function to count number of posts and average per day during academic time and vacation time

analyze_posts <- function(academic_year_filter) {
  # Filtering the data based on the academic_year_filter
  filtered_data <- data_posts %>%
    filter(academic_year == academic_year_filter)
  
  # Count of days for each HEI
  unique_days <- filtered_data %>%
    group_by(id) %>%
    summarise(unique_days = n_distinct(as.Date(created_at)))
  
  # Count of posts for each HEI
  number_posts_boolean <- filtered_data %>%
    group_by(id) %>%
    summarise(count = n())
  
  # Naming the column respecting the time period
  time <- ifelse(academic_year_filter, "academic_time", "vacation_time")
  column_name <- paste0("avg_posts_in_", time)
  
  # Combination of data and calculation of average posts per day
  combined_data <- left_join(unique_days, number_posts_boolean, by = "id")
  combined_data <- combined_data %>%
    mutate(!!column_name := round((count / unique_days), 2))
  
  print(combined_data)
  return(combined_data)
}
data_posts_academic <- analyze_posts(TRUE)
data_posts_vacations <- analyze_posts(FALSE)

Plot for the average number of posts during academic time for each HEI

barplot(data_posts_academic$avg_posts_in_academic_time,
        names.arg = data_posts_academic$id,
        main = "Average Posts during Academic Time",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(data_posts_academic$avg_posts_in_academic_time) + 5),
        las = 2,
        col = "#34495E")

text(x = barplot(data_posts_academic$avg_posts_in_academic_time, plot = FALSE),
     y = data_posts_academic$avg_posts_in_academic_time,
     labels = round(data_posts_academic$avg_posts_in_academic_time, 2),
     pos = 3)

Plot for the average number of posts during vacation time for each HEI

barplot(data_posts_vacations$avg_posts_in_vacation_time,
        names.arg = data_posts_vacations$id,
        main = "Average Posts during Vacation Time",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(data_posts_vacations$avg_posts_in_vacation_time) + 5),
        las = 2,
        col = "#D35400")

text(x = barplot(data_posts_vacations$avg_posts_in_vacation_time, plot = FALSE),
     y = data_posts_vacations$avg_posts_in_vacation_time,
     labels = round(data_posts_vacations$avg_posts_in_vacation_time, 2),
     pos = 3)

Data preparation for dates

# Creating new table that contains a new column for the day of the week
data_posts_days <- data_posts %>%
  mutate(day_of_week = weekdays(created_at))

# Selecting only the id, created_at, and day_of_week columns for the new table
data_posts_days <- data_posts_days %>%
  select(id, created_at, day_of_week)

# Create column hour from created_at
data_posts_days$created_hour <- as.numeric(format(data_posts_days$created_at, "%H"))

print(data_posts_days)
# Grouping by id and day_of_week, then counting the number of posts
number_posts_days <- data_posts_days %>%
  group_by(id, day_of_week) %>%
  summarise(count = n())
`summarise()` has grouped output by 'id'. You can override using the `.groups` argument.
# Grouping by id, day_of_week and day created at, then counting th enumber of tweets
number_posts_per_day <- data_posts_days %>%
    mutate(created_date = as.Date(created_at)) %>% 
    group_by(id, day_of_week, created_date) %>%
    summarize(count = n())
`summarise()` has grouped output by 'id', 'day_of_week'. You can override using the `.groups` argument.
# Finding for each HEI the average count of posts per day
average_number_posts_per_day <- number_posts_per_day %>%
  group_by(id, day_of_week) %>%
  summarise(average_count = round(mean(count), 2))
`summarise()` has grouped output by 'id'. You can override using the `.groups` argument.
print(number_posts_days)

Highest and lowest posts

# Finding the HEI with the lowest count of posts per day
lowest_count <- number_posts_days %>%
  group_by(day_of_week) %>%
  slice_min(order_by = count) %>%
  select(day_of_week, id, count)

# Finding the HEI with the highest count of posts per day
highest_count <- number_posts_days %>%
  group_by(day_of_week) %>%
  slice_max(order_by = count) %>%
  select(day_of_week, id, count)

# Combine the results
high_low_HEI <- bind_rows(lowest_count, highest_count) %>%
  arrange(day_of_week)

print(high_low_HEI)

Plot for the highest and lowest count of posts per day for each day of the week

ggplot(high_low_HEI, aes(x = day_of_week, y = count, fill = id)) +
  geom_bar(stat = "identity", position = "dodge") +
  geom_text(aes(label = count),
            position = position_dodge(width = 0.9),
            vjust = -0.5,
            size = 3) +
  labs(title = "Highest and Lowest Count of Posts per Day for Each Day of the Week",
       x = "Day of the Week", y = "Count") +
  scale_fill_manual(values = rainbow(length(unique(high_low_HEI$id)))) +
  theme_minimal() +
  theme(legend.title = element_blank())

Average of posts

# Finding the HEI with lowest and highest averaged count of posts per day
high_low_average_HEIs <- average_number_posts_per_day %>%
  group_by(day_of_week) %>%
  filter(average_count == max(average_count) | average_count == min(average_count)) %>%
  arrange(day_of_week, ifelse(average_count == min(average_count), average_count, -average_count))

print(high_low_average_HEIs)

Plot for the highest and lowest average count of posts per day for each day of the week

ggplot(high_low_average_HEIs, aes(x = day_of_week, y = average_count, fill = id)) +
  geom_bar(stat = "identity", position = "dodge") +
  geom_text(aes(label = round(average_count, 2)),
            position = position_dodge(width = 0.7),
            vjust = -0.5,
            size = 3) +
  labs(title = "Highest and Lowest Average Count of Posts per Day for Each Day of the Week",
       x = "Day of the Week", y = "Average Count") +
  scale_fill_manual(values = rainbow(length(unique(high_low_HEI$id)))) +
  theme_minimal() +
  theme(legend.title = element_blank())

Favourite hour and day

favourite_day_hei <- number_posts_days %>%
  group_by(id) %>%
  top_n(1, count) %>%
  arrange(id)

print(favourite_day_hei)
# Calculating the number of posts per hour
number_posts_hours <- data_posts_days %>%
  group_by(id, created_hour) %>%
  summarise(count = n()) %>%
  ungroup()
`summarise()` has grouped output by 'id'. You can override using the `.groups` argument.
# Identifying the favorite hour for each HEI
favourite_hour_hei <- number_posts_hours %>%
  group_by(id) %>%
  top_n(1, count) %>%
  arrange(id)

# Adding new columns to categorize the hours and assign numerical value for clustering
favourite_hour_hei <- favourite_hour_hei %>%
  mutate(time_of_day = case_when(
    created_hour >= 5 & created_hour < 12 ~ "Morning",
    created_hour >= 12 & created_hour < 18 ~ "Afternoon",
    TRUE ~ "Night"
  ),
  time_of_day_value = case_when(
    created_hour >= 5 & created_hour < 12 ~ 1,
    created_hour >= 12 & created_hour < 18 ~ 3,
    TRUE ~ 5
  ))

print(favourite_hour_hei)

Heatmaps

Function to plot heatmap for various HEIs

heatmap_maker <- function(target_id){
  # Filtering data for the specific HEI
  target_data <- data_posts_days %>%
    filter(id == target_id)
  
  # Grouping by day of the week and hour, and counting the number of tweets
  tweet_counts <- target_data %>%
    group_by(day_of_week, created_hour) %>%
    summarise(num_posts = n())
  
  # Plotting heatmap
  ggplot(tweet_counts, aes(x = day_of_week, y = created_hour, fill = num_posts)) +
    geom_tile() +
    scale_fill_gradient(low = "white", high = "blue") +
    labs(title = paste("Post Heatmap for", target_id),
         x = "Day of the week",
         y = "Hour of the day")
}

Plot of heatmap for each HEI

heatmap_maker("duke.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("epfl.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("goe.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("harvard.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("leicester.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("manchester.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("mit.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("sb.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("stanford.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("trinity.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("wv.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("yale.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

Hashtags

# Transforming empty strings into NA
data_posts$hashtags[data_posts$hashtags == ""] <- NA

# Table with number of unique hashtags and percentage of usage
hashtags <- data_posts %>%
                group_by(id) %>%
                summarise(count = n(),
                          na = sum(is.na(hashtags)),
                          unique_hashtags = length(unique(hashtags)),
                          hashtag_percentage = round(((count - na) / count * 100), 2))

print(hashtags)

Plot for the count of unique hashtags for each HEI

barplot(hashtags$unique_hashtags,
        names.arg = hashtags$id,
        main = "Unique Hashtags for Each HEI",
        xlab = "HEI",
        ylab = "Count of Unique Hashtags",
        ylim = c(0, max(hashtags$unique_hashtags) + 50),
        las = 2,
        col= "#16A085")

text(x = barplot(hashtags$unique_hashtags, plot = FALSE),
     y = hashtags$unique_hashtags,
     labels = round(hashtags$unique_hashtags, 2),
     pos = 3)

Plot for the usage of hashtag for each HEI

barplot(hashtags$hashtag_percentage,
        names.arg = hashtags$id,
        main = "Hashtags Percentage for Each HEI",
        xlab = "HEI",
        ylab = "Hashtags Percentage",
        ylim = c(0, max(hashtags$hashtag_percentage) + 30),
        las = 2,
        col= "#F1C40F")

text(x = barplot(hashtags$hashtag_percentage, plot = FALSE),
     y = hashtags$hashtag_percentage,
     labels = round(hashtags$hashtag_percentage, 2),
     pos = 3)

URL usage

# Transforming empty strings into NA
data_posts$urls[data_posts$urls == ""] <- NA

# Table with number of post, number of NA and url percentage of usage
url_usage <- data_posts %>%
                group_by(id) %>%
                summarise(count = n(),
                          na = sum(is.na(urls)),
                          url_percentage = round(((count - na) / count * 100), 2))

print(url_usage)

Plot for the usage of hashtag for each HEI

barplot(url_usage$url_percentage,
        names.arg = url_usage$id,
        main = "Urls Percentage for Each HEI",
        xlab = "HEI",
        ylab = "Urls Percentage",
        ylim = c(0, max(url_usage$url_percentage) + 10),
        las = 2,
        col= "#8E44AD")

text(x = barplot(url_usage$url_percentage, plot = FALSE),
     y = url_usage$url_percentage,
     labels = round(url_usage$url_percentage, 2),
     pos = 3)

Text

data_posts_content <- data_posts %>%
            select(id, text)

# Counting number of words
data_posts_content <- data_posts_content %>%
  mutate(num_words = lengths(strsplit(text, "\\s+")))

# Grouping by HEI and calculate average, minimum, and maximum values of number of words
data_posts_content_metrics <- data_posts_content %>%
  group_by(id) %>%
  summarise(average_num_words = mean(num_words),
            min_num_words = min(num_words),
            max_num_words = max(num_words))
print(data_posts_content_metrics)

Plot for the average, maximum and minimum values of words for each HEI

ggplot(data_posts_content_metrics, aes(x = id, y = average_num_words)) +
  geom_point(aes(color = "Average")) +
  geom_errorbar(aes(ymin = min_num_words, ymax = max_num_words, color = "Range"), width = 0.2) +
  scale_color_manual(values = c("Average" = "#1976D2", "Range" = "#EF5350")) +
  labs(title = "Word Count Summary by HEI",
       x = "HEI",
       y = "Number of Words",
       color = "Metric") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust = 1))

Posts Classification

cleanup <- function(docs, spec.words=NULL){
  # lowercase
  docs <- tm_map(docs, content_transformer(tolower))
  # rm numbers
  docs <- tm_map(docs, removeNumbers)
  # rm english common stopWords
  docs <- tm_map(docs, removeWords, stopwords("english"))
  # if stopwords are specified as a character vector
  if(!is.null(spec.words))
    docs <- tm_map(docs, removeWords, spec.words)
  # rm punctuations
  docs <- tm_map(docs, removePunctuation)
  # rm extra white spaces
  docs <- tm_map(docs, stripWhitespace)
  # lemmatizing text
  docs <- tm_map(docs, lemmatize_words)
  
  docs
}
pairing <- function(word_dict, paired_words) {
  if (nrow(word_dict) != length(paired_words)) {
    stop("The number of rows in word_dict and the length of paired_words should be equal.")
  }
  
  word_dict <- tibble::tibble(word = word_dict$value)
  result <- tibble::tibble(word = word_dict$word, category = paired_words)
  return(result)
}
classify_text <- function(text, word_pair) {
  # Tokenizing the text and converting to lowercase
  words <- tolower(unlist(strsplit(text, "\\W+")))
  
  # Finding the frequent words in the text
  freq_words <- words[words %in% word_pair$word]
  
  if (length(freq_words) == 0) {
    return("Unknown")  # Returning Unknown if no frequent words are found
  }
  
  # Getting the corresponding categories
  categories <- word_pair$category[word_pair$word %in% freq_words]
  
  # Sorting categories by frequency in descending order and return the most frequent one
  return(names(sort(table(categories), decreasing = TRUE))[1])
}
corpus_maker <- function(text) {
  texts <- text$text
  vc <- VectorSource(texts)
  corpus <- Corpus(vc)
  
  return(corpus)
}
freq_terms <- function(clean_text, number) {
  dtm <- DocumentTermMatrix(clean_text)
  dtm.tfidf <- weightTfIdf(dtm)
  
  mdf <- as_tibble(as.matrix(dtm.tfidf))
  
  mdf.freq <- mdf %>%
  select(findFreqTerms(dtm, number)) %>%
  summarise_all(sum) %>%
  gather() %>%
  arrange(desc(value))

  mdf.freq$key <- 
    factor(mdf.freq$key,
           levels = mdf.freq$key[order(mdf.freq$value)])
  
  word_dictionary <- as_tibble(mdf.freq$key)
  
  ggplot(mdf.freq, aes(x=key, y=value)) +
    geom_bar(stat="identity") +
    labs(x="terms", y="freq") + coord_flip()
  
  return(word_dictionary)
}

Duke

duke_text <- subset(data_posts_content, id == "duke.csv")
duke_corpus <- corpus_maker(duke_text)

duke_clean_text <-cleanup(duke_corpus, c("new", "will", "change", "north", "can", "first", "year", "carolina", "years", "summer", "meet", "one"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
duke_word_dictionary <- freq_terms(duke_clean_text, 30)

duke_word_category <- c("Image", "Education", "Education", "Image", "Image", "Research", "Research", "Image", "Image", "Engagement", "Research", "Research", "Image", "Education", "Education", "Society", "Education", "Engagement")

duke_word_pair <- pairing(duke_word_dictionary, duke_word_category)

print(duke_word_pair)

# Applying the classify_text function to each text
duke_text <- duke_text %>%
  mutate(category = sapply(text, classify_text, word_pair = duke_word_pair))

print(duke_text)

epfl

epfl_text <- subset(data_posts_content, id == "epfl.csv")
epfl_corpus <- corpus_maker(epfl_text)

epfl_clean_text <- cleanup(epfl_corpus, c("new", "amp", "can", "will", "one", "now", "–", "via", "portrait"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
epfl_word_dictionary <- freq_terms(epfl_clean_text, 30)

epfl_word_category <- c("Image", "Research", "Image", "Research", "Research", "Research", "Education", "Image", "Research", "Research", "Education", "Research", NA, "Education", "Education")

epfl_word_pair <- pairing(epfl_word_dictionary, epfl_word_category)

print(epfl_word_pair)

# Applying the classify_text function to each text
epfl_text <- epfl_text %>%
  mutate(category = sapply(text, classify_text, word_pair = epfl_word_pair))

print(epfl_text)

goe

goe_text <- subset(data_posts_content, id == "goe.csv")
goe_corpus <- corpus_maker(goe_text)

goe_clean_text <- cleanup(goe_corpus, c("can", "new", "amp", "will"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
goe_word_dictionary <- freq_terms(goe_clean_text, 15)

goe_word_category <- c("Research", "Research", "Image", "Education", "Engagement", "Image", "Engagement", "Education")

goe_word_pair <- pairing(goe_word_dictionary, goe_word_category)

print(goe_word_pair)

# Applying the classify_text function to each text
goe_text <- goe_text %>%
  mutate(category = sapply(text, classify_text, word_pair = goe_word_pair))

print(goe_text)

harvard

harvard_text <- subset(data_posts_content, id == "harvard.csv")
harvard_corpus <- corpus_maker(harvard_text)

harvard_clean_text <- cleanup(harvard_corpus, c("new", "can", "will", "summer", "year", "first", "may", "-", "recent", "years", "one", "said", "time", "many", "world", "change"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
harvard_word_dictionary <- freq_terms(harvard_clean_text, 50)

harvard_word_category <- c("Image", "Education", "Research", "Image", "Research", "Research", "Engagement", "Education", "Society", "Education", "Society", "Research", "Education", "Image", "Research", NA, "Research", "Education", "Education", NA, "Research", "Society", "Society", "Research", "Education")

harvard_word_pair <- pairing(harvard_word_dictionary, harvard_word_category)

print(harvard_word_pair)

harvard_text <- harvard_text %>%
  mutate(category = sapply(text, classify_text, word_pair = harvard_word_pair))

print(harvard_text)

leicester

leicester_text <- subset(data_posts_content, id == "leicester.csv")
leicester_corpus <- corpus_maker(leicester_text)

leicester_clean_text <- cleanup(leicester_corpus, c("👉", "day", "new", "clear", "year", "can", "will", "time", "space", "one", "first"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
leicester_word_dictionary <- freq_terms(leicester_clean_text, 60)

leicester_word_category <- c("Image", NA, "Image", "Engagement", "Society", "Education", "Education", "Research", "Engagement", "Engagement", "Image", "Engagement", "Education", "Engamement", "Engagement", "Education", "Image", "Engagement", "Education", NA, "Image")

leicester_word_pair <- pairing(leicester_word_dictionary, leicester_word_category)

print(leicester_word_pair)

leicester_text <- leicester_text %>%
  mutate(category = sapply(text, classify_text, word_pair = leicester_word_pair))

print(leicester_text)

manchester

manchester_text <- subset(data_posts_content, id == "manchester.csv")
manchester_corpus <- corpus_maker(manchester_text)

manchester_clean_text <- cleanup(manchester_corpus, c("👇", "can", "just", "will", "get", "well", "one", "help", "now", "new", "read", "congratulations"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
manchester_word_dictionary <- freq_terms(manchester_clean_text, 50)

manchester_word_category <- c(NA, "Image", "Engagement", "Image", "Education", "Education", "Engagement", "Image", "Education", "Society", "Research", "Research", "Education", "Education", "Education", "Engagement", "Research", "Education", "Research")

manchester_word_pair <- pairing(manchester_word_dictionary, manchester_word_category)

print(manchester_word_pair)

manchester_text <- manchester_text %>%
  mutate(category = sapply(text, classify_text, word_pair = manchester_word_pair))

print(manchester_text)

mit

mit_text <- subset(data_posts_content, id == "mit.csv")
mit_corpus <- corpus_maker(mit_text)

mit_clean_text <- cleanup(mit_corpus, c("new", "can", "says", "“", "’", "—", "", "may", "first", "will", "using", "way", "one", "science"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
mit_word_dictionary <- freq_terms(mit_clean_text, 40)

mit_word_category <- c("Image", "Research", NA, "Image", NA, "Education", "Research", NA, "Education", "Research", "Research", "Education", "Research", "Research", "Image", "Education", "Education", "Research", "Research", "Research", "Research", "Research", "Research", "Research", "Research", "Society")

mit_word_pair <- pairing(mit_word_dictionary, mit_word_category)

print(mit_word_pair)

mit_text <- mit_text %>%
  mutate(category = sapply(text, classify_text, word_pair = mit_word_pair))

print(mit_text)

sb

sb_text <- subset(data_posts_content, id == "sb.csv")
sb_corpus <- corpus_maker(sb_text)

sb_clean_text <- cleanup(sb_corpus, c("new", "—", "will", "week", "can", "amp", "via", "now", "future", "first"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
sb_word_dictionary <- freq_terms(sb_clean_text, 30)

sb_word_category <- c("Image", "Research", "Image", "Image", "Image", NA, "Research", "Image", "Image", "Research", "Education", "Image", "Education", "Image", "Education", "Image", "Society")

sb_word_pair <- pairing(sb_word_dictionary, sb_word_category)

print(sb_word_pair)

sb_text <- sb_text %>%
  mutate(category = sapply(text, classify_text, word_pair = sb_word_pair))

print(sb_text)

stanford

stanford_text <- subset(data_posts_content, id == "stanford.csv")
stanford_corpus <- corpus_maker(stanford_text)

stanford_clean_text <- cleanup(stanford_corpus, c("new", "will"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
stanford_word_dictionary <- freq_terms(stanford_clean_text, 25)

stanford_word_category <- c("Image", "Image", "Research", "Education", NA)

stanford_word_pair <- pairing(stanford_word_dictionary, stanford_word_category)

print(stanford_word_pair)

stanford_text <- stanford_text %>%
  mutate(category = sapply(text, classify_text, word_pair = stanford_word_pair))

print(stanford_text)

trinity

trinity_text <- subset(data_posts_content, id == "trinity.csv")
trinity_corpus <- corpus_maker(trinity_text)

trinity_clean_text <- cleanup(trinity_corpus, c("amp", "read", "new", "can", "will", "week", "work", "great", "day", "visit", "irish", "first", "led", "congratulations"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
trinity_word_dictionary <- freq_terms(trinity_clean_text, 40)
Warning: empty document(s): 754
trinity_word_category <- c("Image", "Education", "Research", "Image", "Research", "Image", "Research", "Society", "Society", "Education", "Engagement", "Engagement", "Education", "Image", "Education", "Image", "Research")

trinity_word_pair <- pairing(trinity_word_dictionary, trinity_word_category)

print(trinity_word_pair)

trinity_text <- trinity_text %>%
  mutate(category = sapply(text, classify_text, word_pair = trinity_word_pair))

print(trinity_text)

wv

wv_text <- subset(data_posts_content, id == "wv.csv")
wv_corpus <- corpus_maker(wv_text)

wv_clean_text <- cleanup(wv_corpus, c("💛💙", "🙌", "👉", "happy", "day", "see", "week", "great", "can", "will", "well", "know", "new", "now", "get", "just"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
wv_word_dictionary <- freq_terms(wv_clean_text, 40)

wv_word_category <- c(NA, "Image", "Image", NA, "Image", "Image", "Engagement", "Education", "Education", "Education", NA)

wv_word_pair <- pairing(wv_word_dictionary, wv_word_category)

print(wv_word_pair)

wv_text <- wv_text %>%
  mutate(category = sapply(text, classify_text, word_pair = wv_word_pair))

print(wv_text)

yale

yale_text <- subset(data_posts_content, id == "yale.csv")
yale_corpus <- corpus_maker(yale_text)

yale_clean_text <- cleanup(yale_corpus, c("new", "—", "will", "can", "'", "first", "work", "read", "help", "year"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
yale_word_dictionary <- freq_terms(yale_clean_text, 70)

yale_word_category <- c("Research", "Education", "Education", NA, "Research", "Research", "Research", "Research", "Education", "Image", "Research", "Education", "Image", "Research", NA, "Research", "Image", "Image", "Research", "Image", "Research", "Society", "Research", "Society", "Society")

yale_word_pair <- pairing(yale_word_dictionary, yale_word_category)

print(yale_word_pair)

yale_text <- yale_text %>%
  mutate(category = sapply(text, classify_text, word_pair = yale_word_pair))

print(yale_text)

Sentiment Analysis

Function for emotions

emotions_maker <- function(texts, hei_name){
  text_emotion <- get_nrc_sentiment(texts)
  
  # Proportions for text values 
  emotion_proportions <- colSums(prop.table(text_emotion[, 1:8]))
  
  # Visualization of percentages on each emotion found on posts
  barplot <- barplot(
    sort(colSums(prop.table(text_emotion[, 1:8]))), 
    horiz = TRUE, 
    cex.names = 0.7, 
    las = 1, 
    main = sprintf("Emotions found in %s's texts", hei_name), 
    xlab="Percentage",
    xlim = c(0, max(emotion_proportions) * 1.2)
  )
  
  text(
    x = sort(emotion_proportions),
    y = barplot,                    
    labels = sprintf("%.2f%%", 100 * sort(emotion_proportions)),
    pos = 4,                               
    cex = 0.7                              
  )
}

Duke

duke_text$sentiment <- round(get_sentiment(duke_text$text), 2)

print(duke_text)
emotions_maker(duke_text$text, "duke")

epfl

epfl_text$sentiment <- round(get_sentiment(epfl_text$text), 2)

print(epfl_text)
emotions_maker(epfl_text$text, "epfl")

goe

goe_text$sentiment <- round(get_sentiment(goe_text$text), 2)

print(goe_text)
emotions_maker(goe_text$text, "goe")

harvard

harvard_text$sentiment <- round(get_sentiment(harvard_text$text), 2)

print(harvard_text)
emotions_maker(harvard_text$text, "harvard")

leicester

leicester_text$sentiment <- round(get_sentiment(leicester_text$text), 2)

print(leicester_text)
emotions_maker(leicester_text$text, "leicester")

manchester

manchester_text$sentiment <- round(get_sentiment(manchester_text$text), 2)

print(manchester_text)
emotions_maker(manchester_text$text, "manchester")

mit

mit_text$sentiment <- round(get_sentiment(mit_text$text), 2)

print(mit_text)
emotions_maker(mit_text$text, "mit")

sb

sb_text$sentiment <- round(get_sentiment(sb_text$text), 2)

print(sb_text)
emotions_maker(sb_text$text, "sb")

stanford

stanford_text$sentiment <- round(get_sentiment(stanford_text$text), 2)

print(stanford_text)
emotions_maker(stanford_text$text, "stanford")

trinity

trinity_text$sentiment <- round(get_sentiment(trinity_text$text), 2)

print(trinity_text)
emotions_maker(trinity_text$text, "trinity")

wv

wv_text$sentiment <- round(get_sentiment(wv_text$text), 2)

print(wv_text)
emotions_maker(wv_text$text, "wv")

yale

yale_text$sentiment <- round(get_sentiment(yale_text$text), 2)

print(yale_text)
emotions_maker(yale_text$text, "yale")

Creation of average sentiment for each HEI

# Combine all the sentiment scores into one data frame
all_sentiments <- rbind(
    duke_text[, c("id", "sentiment")],
    epfl_text[, c("id", "sentiment")],
    goe_text[, c("id", "sentiment")],
    harvard_text[, c("id", "sentiment")],
    leicester_text[, c("id", "sentiment")],
    manchester_text[, c("id", "sentiment")],
    mit_text[, c("id", "sentiment")],
    sb_text[, c("id", "sentiment")],
    stanford_text[, c("id", "sentiment")],
    trinity_text[, c("id", "sentiment")],
    wv_text[, c("id", "sentiment")],
    yale_text[, c("id", "sentiment")]
)

# Calculate the average sentiment for each HEI
hei_average_sentiments <- aggregate(sentiment ~ id, data = all_sentiments, FUN = mean) %>%
    rename(average_sentiment = sentiment)

# Print the table
print(hei_average_sentiments)

Clusters

Function for cosine matrix

cosine_matrix_maker <- function(table){
  numerical_columns <- table %>%
    select(-id)
  
  # Normalizing columns
  normalized_columns <- as.data.frame(scale(numerical_columns))
  
  # Transposing the data to compute similarity between columns
  transposed_data <- t(normalized_columns)
  
  # Computing the cosine similarity matrix
  similarity_matrix <- as.matrix(proxy::dist(transposed_data, method = "cosine"))
  
  # Converting distance to similarity
  similarity_matrix <- 1 - similarity_matrix
  
  # Print the similarity matrix
  print(similarity_matrix)
}  

Creating table for cluster algorithms

# Joining attribute percentage_tweets (percentage of tweets out of all posts) and percentage_replies (percentage of replies out of all posts) from number_posts also adding unique_hashtags (number of unique hashtags) and hashtag_percentage (percentage of posts that contain a hashtag) from hashtags, per HEI
cluster_table <- merge(select(hashtags, id, unique_hashtags, hashtag_percentage), select(data_ratio, id, percentage_tweets, percentage_replies), by = "id", all=TRUE)

# Joining attribute avg_posts_per_days (average of posts per day) from posts_per_day per HEI
cluster_table <- merge(cluster_table, select(posts_per_day, id, avg_posts_per_days), by = "id", all=TRUE)

# Joining attribute avg_posts_per_weeks (average of posts per week) from posts_per_week per HEI
cluster_table <- merge(cluster_table, select(posts_per_week, id, avg_posts_per_weeks), by = "id", all=TRUE)

# Joining attribute avg_posts_in_academic_time (average of posts during academic time)  from data_posts_academic per HEI
cluster_table <- merge(cluster_table, select(data_posts_academic, id, avg_posts_in_academic_time), by = "id", all=TRUE)

# Joining attribute avg_posts_in_vacation_time (average of posts during vacation time) from data_posts_vacations per HEI
cluster_table <- merge(cluster_table, select(data_posts_vacations, id, avg_posts_in_vacation_time), by = "id", all=TRUE)

# Joining attribute time_of_day_value (numerical value referring to time of day where every HEI made more posts) from favourite_hour_hei per HEI
cluster_table <- merge(cluster_table, select(favourite_hour_hei, id, time_of_day_value), by = "id", all=TRUE)

# Joining attribute url_percentage (percentage of posts that contain an url) from url_usage per HEI
cluster_table <- merge(cluster_table, select(url_usage, id, url_percentage), by = "id", all=TRUE)

# Joining attribute average_num_words (average number of words in the posts) from data_posts_content_metrics per HEI
cluster_table <- merge(cluster_table, select(data_posts_content_metrics, id, average_num_words), by = "id", all=TRUE)

# Joining attribute average_sentiment (average sentiment of posts) from hei_average_sentiments per HEI
cluster_table <- merge(cluster_table, select(hei_average_sentiments, id, average_sentiment), by = "id", all=TRUE)

print(cluster_table)
cosine_matrix_maker(cluster_table)
                           unique_hashtags hashtag_percentage percentage_tweets percentage_replies avg_posts_per_days avg_posts_per_weeks
unique_hashtags                 1.00000000        0.714423178        0.16552198        -0.15047782         0.09423180          0.09460315
hashtag_percentage              0.71442318        1.000000000        0.13470812        -0.34565768        -0.01048208         -0.01015282
percentage_tweets               0.16552198        0.134708118        1.00000000        -0.73554060         0.30916774          0.30775916
percentage_replies             -0.15047782       -0.345657681       -0.73554060         1.00000000         0.09388665          0.09480553
avg_posts_per_days              0.09423180       -0.010482084        0.30916774         0.09388665         1.00000000          0.99999797
avg_posts_per_weeks             0.09460315       -0.010152816        0.30775916         0.09480553         0.99999797          1.00000000
avg_posts_in_academic_time      0.12447464       -0.033819668        0.32273419         0.06234258         0.97881722          0.97877880
avg_posts_in_vacation_time      0.09368691        0.004006656        0.16186441         0.20825601         0.97718481          0.97729844
time_of_day_value              -0.22814041       -0.175920221        0.30874190        -0.21552038        -0.05609489         -0.05648404
url_percentage                  0.09757172        0.395409742        0.52625841        -0.84342954         0.06010475          0.05920162
average_num_words               0.46150077        0.713419420        0.06039969        -0.48123350        -0.36976005         -0.36993229
average_sentiment               0.59299320        0.645971142       -0.31588590        -0.04093998        -0.49647197         -0.49605792
                           avg_posts_in_academic_time avg_posts_in_vacation_time time_of_day_value url_percentage average_num_words average_sentiment
unique_hashtags                            0.12447464                0.093686911       -0.22814041     0.09757172        0.46150077        0.59299320
hashtag_percentage                        -0.03381967                0.004006656       -0.17592022     0.39540974        0.71341942        0.64597114
percentage_tweets                          0.32273419                0.161864406        0.30874190     0.52625841        0.06039969       -0.31588590
percentage_replies                         0.06234258                0.208256011       -0.21552038    -0.84342954       -0.48123350       -0.04093998
avg_posts_per_days                         0.97881722                0.977184814       -0.05609489     0.06010475       -0.36976005       -0.49647197
avg_posts_per_weeks                        0.97877880                0.977298444       -0.05648404     0.05920162       -0.36993229       -0.49605792
avg_posts_in_academic_time                 1.00000000                0.953796153       -0.14950356     0.07310820       -0.37164829       -0.49875096
avg_posts_in_vacation_time                 0.95379615                1.000000000       -0.07107087     0.02168275       -0.32110128       -0.41192255
time_of_day_value                         -0.14950356               -0.071070873        1.00000000    -0.01251319       -0.10921047       -0.29227459
url_percentage                             0.07310820                0.021682748       -0.01251319     1.00000000        0.63011343        0.17847188
average_num_words                         -0.37164829               -0.321101278       -0.10921047     0.63011343        1.00000000        0.84649940
average_sentiment                         -0.49875096               -0.411922553       -0.29227459     0.17847188        0.84649940        1.00000000

Based on this cosine matrix we decided to remove average_num_words (high similarity to unique_hashtags, hashtag_percentage and url_percentage), unique_hashtags (high similarity to hashtag_percentage and the distances of hashtag_percentage to the other columns is lower) and avg_posts_per_weeks, avg_posts_in_academic_time, avg_posts_in_vacation_time are also all removed (due to their similarity to avg_posts_per_days)

cluster_table <- cluster_table %>%
  select(-average_num_words, -unique_hashtags, -avg_posts_per_weeks, -avg_posts_in_academic_time, -avg_posts_in_vacation_time)

cosine_matrix_maker(cluster_table)
                   hashtag_percentage percentage_tweets percentage_replies avg_posts_per_days time_of_day_value url_percentage average_sentiment
hashtag_percentage         1.00000000         0.1347081        -0.34565768        -0.01048208       -0.17592022     0.39540974        0.64597114
percentage_tweets          0.13470812         1.0000000        -0.73554060         0.30916774        0.30874190     0.52625841       -0.31588590
percentage_replies        -0.34565768        -0.7355406         1.00000000         0.09388665       -0.21552038    -0.84342954       -0.04093998
avg_posts_per_days        -0.01048208         0.3091677         0.09388665         1.00000000       -0.05609489     0.06010475       -0.49647197
time_of_day_value         -0.17592022         0.3087419        -0.21552038        -0.05609489        1.00000000    -0.01251319       -0.29227459
url_percentage             0.39540974         0.5262584        -0.84342954         0.06010475       -0.01251319     1.00000000        0.17847188
average_sentiment          0.64597114        -0.3158859        -0.04093998        -0.49647197       -0.29227459     0.17847188        1.00000000

Function for cluster method

cluster_maker <- function(num_clusters, table){
  # Excluding id column for clustering
  cluster_data <- select(table, -id)
    
  # Scaling the data for kmeans method
  scaled_data <- scale(cluster_data)
  
  kmeans_model <- kmeans(scaled_data, centers = num_clusters, nstart = 10)

  # Extract cluster assignments
  cluster_assignments <- kmeans_model$cluster
  
  # Create a data frame combining original data with cluster assignments
  clustered_data <- cbind(cluster_table$id, cluster_data, cluster = cluster_assignments)
  
  clustered_data <- clustered_data[, c("cluster_table$id", "cluster")]
  
  print(clustered_data)
}

Function to discover best number of clusters

elbow_maker <- function(table){
  cluster_data <- select(table, -id)
  scaled_data <- scale(cluster_data)
  
  wss <- vector()
  range <- 1:10
  
  for (k in range) {
    kmeans_model <- kmeans(scaled_data, centers = k, nstart = 10)
    wss[k] <- kmeans_model$tot.withinss
  }
  
  elbow_df <- data.frame(k = range, WSS = wss)
  ggplot(elbow_df, aes(x = k, y = WSS)) +
    geom_line() +
    geom_point() +
    labs(x = "Number of Clusters", y = "Within-Cluster Sum of Squares (WCSS)",
         title = "Elbow Method for Optimal k") +
    theme_minimal()
}

Plot of Elbow Method and selection of best number of cluster to view how HEIs are grouped

elbow_maker(cluster_table)

cluster_maker(4, cluster_table)
LS0tCnRpdGxlOiAiUHJvamVjdCIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQojIyMgRm9yIFIgYmVnaW5uZXJzCk5ldyBjaHVuayAqQ3RybCtBbHQrSSoKCkV4ZWN1dGUgY2h1bmsgKkN0cmwrU2hpZnQrRW50ZXIqCgpFeGVjdXRlIGFsbCBjaHVua3MgKkN0cmwrQWx0K1IqCgpIVE1MIHByZXZpZXcgKkN0cmwrU2hpZnQrSyoKCiMgTGlicmFyeSBwcmVwYXJhdGlvbnMKCmBgYHtyfQpsaWJyYXJ5KHJlYWRyKQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KHJlc2hhcGUyKQpsaWJyYXJ5KHN0YXRzKQpsaWJyYXJ5KHRtKQpsaWJyYXJ5KHRleHQydmVjKQpsaWJyYXJ5KHRleHRzdGVtKQpsaWJyYXJ5KHN5dXpoZXQpCmBgYAoKIyBEYXRhIEltcG9ydAoKYGBge3J9CmRhdGFfcG9zdHMgPC0gcmVhZC5jc3YoIn4vNHllYXIvMnNlbWVzdGVyL2R0SUkvQ1NWcy9IRUlzLmNzdiIsCiAgICAgICAgICAgICAgICAgY29sQ2xhc3NlcyA9IGModHdlZXRfaWQgPSAiY2hhcmFjdGVyIikpCgojIE1vZGlmeWluZyBjcmVhdGVkX2F0IHR5cGUgc28gdGhhdCBhdHRyaWJ1dGUgY2FuIGJlIHVzZWQgbW9yZSBlYXNpbHkgCmRhdGFfcG9zdHMkY3JlYXRlZF9hdCA8LSBhcy5QT1NJWGN0KGRhdGFfcG9zdHMkY3JlYXRlZF9hdCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZm9ybWF0PSAiJVktJW0tJWRUJUg6JU06JVMiLCB0ej0iVVRDIikKCiNWaWV3KGRhdGEpCnN1bW1hcnkoZGF0YV9wb3N0cykKYGBgCgojIEluaXRpYWwgRGF0YSBQcmVwYXJhdGlvbgoKYGBge3J9CiMgQ291bnQgb2YgaG93IG1hbnkgZW50cmllcyBlYWNoIEhFSSBoYXMKbnVtYmVyX3Bvc3RzIDwtIGRhdGFfcG9zdHMgJT4lCiAgICAgICAgICAgICAgZ3JvdXBfYnkoaWQpICU+JSBzdW1tYXJpc2UoY291bnQgPSBuKCkpCgpudW1iZXJfcG9zdHMKYGBgCgojIFNpbmNlIGNvbXBsdXRlbnNlIG9ubHkgaGFzIDEgZW50cnkgd2UgY2FuJ3QgbGVhcm4gYW55dGhpbmcgZnJvbSBpdCwgc28gd2UgcmVtb3ZlZCBpdAoKYGBge3J9CmRhdGFfcG9zdHMgPC0gZGF0YV9wb3N0c1tkYXRhX3Bvc3RzJGlkICE9ICJjb21wbHV0ZW5zZS5jc3YiLCBdCmBgYAoKIyBWaXN1YWxpemF0aW9uIG9mIGFsbCBwb3N0cywganVzdCB0d2VldHMgYW5kIGp1c3QgcmVwbGllcwoKYGBge3J9Cm51bWJlcl9wb3N0cyA8LSBkYXRhX3Bvc3RzICU+JQogICAgICAgICAgICAgIGdyb3VwX2J5KGlkKSAlPiUgc3VtbWFyaXNlKHBvc3RzID0gbigpKQoKbnVtYmVyX3R3ZWV0cyA8LSBkYXRhX3Bvc3RzW2RhdGFfcG9zdHMkdHlwZSA9PSAiVHdlZXQiLCBdICU+JQogICAgICAgICAgICAgIGdyb3VwX2J5KGlkKSAlPiUgc3VtbWFyaXNlKHR3ZWV0cyA9IG4oKSkKCm51bWJlcl9yZXBsaWVzIDwtIGRhdGFfcG9zdHNbZGF0YV9wb3N0cyR0eXBlID09ICJSZXBseSIsIF0gJT4lCiAgICAgICAgICAgICAgZ3JvdXBfYnkoaWQpICU+JSBzdW1tYXJpc2UocmVwbGllcyA9IG4oKSkKCnByaW50KG51bWJlcl9wb3N0cykKcHJpbnQobnVtYmVyX3R3ZWV0cykKcHJpbnQobnVtYmVyX3JlcGxpZXMpCmBgYAoKIyBDYWxjdWxhdGluZyB0aGUgcGVyY2VudGFnZSBvZiB0d2VldHMgYW5kIHJlcGxpZXMgYmFzZWQgb24gYWxsIHBvc3RzCgpgYGB7cn0KIyBNZXJnaW5nIHRoZSBjb3VudHMgb2YgdHdlZXRzIGFuZCByZXBsaWVzIHdpdGggdGhlIGNvdW50IG9mIHBvc3RzCmRhdGFfcmF0aW8gPC0gbWVyZ2UobnVtYmVyX3Bvc3RzLCBudW1iZXJfdHdlZXRzLCBieSA9ICJpZCIsIGFsbCA9IFRSVUUpCmRhdGFfcmF0aW8gPC0gbWVyZ2UoZGF0YV9yYXRpbywgbnVtYmVyX3JlcGxpZXMsIGJ5ID0gImlkIiwgYWxsID0gVFJVRSkKCgpkYXRhX3JhdGlvJHBlcmNlbnRhZ2VfdHdlZXRzIDwtIHJvdW5kKCgoZGF0YV9yYXRpbyR0d2VldHMgLyBkYXRhX3JhdGlvJHBvc3RzKSAqIDEwMCksIDIpCmRhdGFfcmF0aW8kcGVyY2VudGFnZV9yZXBsaWVzIDwtIHJvdW5kKCgoZGF0YV9yYXRpbyRyZXBsaWVzIC8gZGF0YV9yYXRpbyRwb3N0cykgKiAxMDApLCAyKQoKZGF0YV9yYXRpbyA8LSBkYXRhX3JhdGlvWywgYygiaWQiLCAicGVyY2VudGFnZV90d2VldHMiLCAicGVyY2VudGFnZV9yZXBsaWVzIildCgpkYXRhX3JhdGlvJHBlcmNlbnRhZ2VfcmVwbGllc1tpcy5uYShkYXRhX3JhdGlvJHBlcmNlbnRhZ2VfcmVwbGllcyldIDwtIDAKCnByaW50KGRhdGFfcmF0aW8pCmBgYAoKIyBOQSByZW1vdmFsCgojIEZ1bmN0aW9uIHRvIHZpc3VhbGl6ZSB0aGUgbnVtYmVyIG9mIE5BcyBpbiBhbGwgY29sdW1ucwoKYGBge3J9Cm5hX2NvdW50IDwtIGZ1bmN0aW9uKCl7CiAgIyBDb3VudGluZyB0aGUgbnVtYmVyIG9mIE5BIHZhbHVlcyBmb3IgZWFjaCBjb2x1bW4KICBuYV9jb3VudCA8LSBjb2xTdW1zKGlzLm5hKGRhdGFfcG9zdHMpKQogIAogICMgQ3JlYXRpbmcgYSBuZXcgZGF0YSBmcmFtZSB3aXRoIHRoZSBOQSBjb3VudHMKICBuYV9jb3VudHNfdGFibGUgPC0gZGF0YS5mcmFtZShDb2x1bW4gPSBuYW1lcyhuYV9jb3VudCksIE5BX0NvdW50ID0gbmFfY291bnQpCiAgCiAgcHJpbnQobmFfY291bnRzX3RhYmxlKQp9CmBgYAoKIyBDYWxjdWxhdGlvbnMgb2YgdmlldywgZmF2b3VyaXRlLCByZXR3ZWV0IGFuZCByZXBseSBwZXJjZW50aWxlcyBhbmQgdmlzdWFsaXphdGlvbiBvZiBOQXMgaW4gYWxsIGNvbHVtbnMKCmBgYHtyfQpkYXRhX3Bvc3RzIDwtIGRhdGFfcG9zdHMgJT4lCiAgZ3JvdXBfYnkoaWQpICU+JQogIG11dGF0ZSh2aWV3X3BlcmNlbnRpbGUgPSBudGlsZSh2aWV3X2NvdW50LCAxMDApLAogICAgICAgICBmYXZvcml0ZV9wZXJjZW50aWxlID0gbnRpbGUoZmF2b3JpdGVfY291bnQsIDEwMCksCiAgICAgICAgIHJldHdlZXRfcGVyY2VudGlsZSA9IG50aWxlKHJldHdlZXRfY291bnQsIDEwMCksCiAgICAgICAgIHJlcGx5X3BlcmNlbnRpbGUgPSBudGlsZShyZXBseV9jb3VudCwgMTAwKSkgJT4lCiAgcm93d2lzZSgpICU+JQogIG11dGF0ZShhdmdfcGVyY2VudGlsZSA9IHJvdW5kKG1lYW4oYyh2aWV3X3BlcmNlbnRpbGUsIGZhdm9yaXRlX3BlcmNlbnRpbGUsIHJldHdlZXRfcGVyY2VudGlsZSwgcmVwbHlfcGVyY2VudGlsZSksIG5hLnJtID0gVFJVRSksIDIpKQoKbmFfY291bnQoKQoKZGF0YV9wZXJjZW50aWxlIDwtIGRhdGFfcG9zdHNbLCBjKCJpZCIsICJ2aWV3X3BlcmNlbnRpbGUiLCAiZmF2b3JpdGVfcGVyY2VudGlsZSIsICJyZXR3ZWV0X3BlcmNlbnRpbGUiLCAicmVwbHlfcGVyY2VudGlsZSIsICJhdmdfcGVyY2VudGlsZSIpXQoKcHJpbnQoZGF0YV9wZXJjZW50aWxlKQpgYGAKCiMgQ2FsY3VsYXRpb24gb2YgdGhlIG1heGltdW0gbnVtYmVyIG9mIHZpZXdzIGZvciBlYWNoIEhFSQoKYGBge3J9Cm1heF92aWV3X2NvdW50cyA8LSB0YXBwbHkoZGF0YV9wb3N0cyR2aWV3X2NvdW50LCBkYXRhX3Bvc3RzJGlkLCBtYXgsIG5hLnJtID0gVFJVRSkKCnByaW50KG1heF92aWV3X2NvdW50cykKYGBgCgojIFJlbW92YWwgb2YgTkFzCgpgYGB7cn0KIyBGcm9tIHZpZXcgY291bnQKZGF0YV9wb3N0cyR2aWV3X2NvdW50IDwtIGlmZWxzZSgKICBpcy5uYShkYXRhX3Bvc3RzJHZpZXdfY291bnQpLAogIHJvdW5kKG1heF92aWV3X2NvdW50c1tkYXRhX3Bvc3RzJGlkXSAqIChkYXRhX3Bvc3RzJGF2Z19wZXJjZW50aWxlIC8gMTAwKSksCiAgZGF0YV9wb3N0cyR2aWV3X2NvdW50KQoKIyBGcm9tIHZpZXcgcGVyY2VudGlsZQpkYXRhX3Bvc3RzJHZpZXdfcGVyY2VudGlsZSA8LSBpZmVsc2UoCiAgaXMubmEoZGF0YV9wb3N0cyR2aWV3X3BlcmNlbnRpbGUpLAogIGRhdGFfcG9zdHMkYXZnX3BlcmNlbnRpbGUsCiAgZGF0YV9wb3N0cyR2aWV3X3BlcmNlbnRpbGUpCmBgYAoKIyBWaXN1YWxpemF0aW9uIG9mIE5BcyBpbiBhbGwgY29sdW1ucwoKYGBge3J9Cm5hX2NvdW50KCkKYGBgCgojIEZ1bmN0aW9uIHRvIGNhbGN1bGF0ZSBhdmVyYWdlIHBvc3RzCgpgYGB7cn0KYXZlcmFnZV9wb3N0cyA8LSBmdW5jdGlvbih0aW1lZnJhbWUpewogICMgQ2FsY3VsYXRpb24gb2YgdGhlIHRpbWVmcmFtZSBiZXR3ZWVuIGVhcmxpZXN0IGFuZCBsYXRlc3QgcG9zdCBmb3IgZWFjaCBIRUkKICBkYXRlX3JhbmdlIDwtIGRhdGFfcG9zdHMgJT4lCiAgICBncm91cF9ieShpZCkgJT4lCiAgICBzdW1tYXJpc2UobWluX2RhdGUgPSBtaW4oY3JlYXRlZF9hdCksCiAgICAgICAgICAgICAgbWF4X2RhdGUgPSBtYXgoY3JlYXRlZF9hdCkpICU+JQogICAgbXV0YXRlKG51bV9kYXlzID0gYXMubnVtZXJpYyhkaWZmdGltZShtYXhfZGF0ZSwgbWluX2RhdGUsIHVuaXRzID0gdGltZWZyYW1lKSkpCiAgCiAgIyBOYW1pbmcgdGhlIGNvbHVtbiByZXNwZWN0aW5nIHRoZSB0aW1lZnJhbWUKICBjb2x1bW5fbmFtZSA8LSBwYXN0ZTAoImF2Z19wb3N0c19wZXJfIiwgdGltZWZyYW1lKQogIAogICMgQ2FsY3VsYXRpb24gb2YgdGhlIG51bWJlciBvZiBwb3N0cyBwZXIgZGF5IGZvciBlYWNoIEhFSQogIHBvc3RzX3Blcl90aW1lZnJhbWUgPC0gbnVtYmVyX3Bvc3RzICU+JQogICAgbGVmdF9qb2luKGRhdGVfcmFuZ2UsIGJ5ID0gImlkIikgJT4lCiAgICBtdXRhdGUoISFjb2x1bW5fbmFtZSA6PSByb3VuZCgocG9zdHMgLyBudW1fZGF5cyksIDIpKQogIAogIHByaW50KHBvc3RzX3Blcl90aW1lZnJhbWUpCiAgcmV0dXJuKHBvc3RzX3Blcl90aW1lZnJhbWUpCn0KYGBgCgpgYGB7cn0KcG9zdHNfcGVyX2RheSA8LSBhdmVyYWdlX3Bvc3RzKCJkYXlzIikKcG9zdHNfcGVyX3dlZWsgPC0gYXZlcmFnZV9wb3N0cygid2Vla3MiKQpgYGAKCiMgUGxvdCBmb3IgdGhlIGF2ZXJhZ2UgbnVtYmVyIG9mIHBvc3RzIHBlciBkYXkgZm9yIGVhY2ggSEVJCgpgYGB7cn0KYmFycGxvdChwb3N0c19wZXJfZGF5JGF2Z19wb3N0c19wZXJfZGF5cywKICAgICAgICBuYW1lcy5hcmcgPSBwb3N0c19wZXJfZGF5JGlkLAogICAgICAgIG1haW4gPSAiQXZlcmFnZSBQb3N0cyBwZXIgRGF5IiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJBdmVyYWdlIE51bWJlciBvZiBQb3N0cyIsCiAgICAgICAgeWxpbSA9IGMoMCwgbWF4KHBvc3RzX3Blcl9kYXkkYXZnX3Bvc3RzX3Blcl9kYXlzKSArIDEpLAogICAgICAgIGxhcyA9IDIsCiAgICAgICAgY29sID0gIiMzNDk4REIiKQoKIyBBZGRpbmcgdGV4dCBsYWJlbHMgb3ZlciBlYWNoIGJhciBhbmQgYWxpZ25pbmcgaXQgd2l0aCB0aGUgY2VudGVyIG9mIGVhY2ggYmFyIAp0ZXh0KHggPSBiYXJwbG90KHBvc3RzX3Blcl9kYXkkYXZnX3Bvc3RzX3Blcl9kYXlzLCBwbG90ID0gRkFMU0UpLAogICAgIHkgPSBwb3N0c19wZXJfZGF5JGF2Z19wb3N0c19wZXJfZGF5cywKICAgICBsYWJlbHMgPSByb3VuZChwb3N0c19wZXJfZGF5JGF2Z19wb3N0c19wZXJfZGF5cywgMiksCiAgICAgcG9zID0gMykKYGBgCgojIFBsb3QgZm9yIHRoZSBhdmVyYWdlIG51bWJlciBvZiBwb3N0cyBwZXIgd2VlayBmb3IgZWFjaCBIRUkKCmBgYHtyfQpiYXJwbG90KHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MsCiAgICAgICAgbmFtZXMuYXJnID0gcG9zdHNfcGVyX3dlZWskaWQsCiAgICAgICAgbWFpbiA9ICJBdmVyYWdlIFBvc3RzIHBlciBXZWVrIiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJBdmVyYWdlIE51bWJlciBvZiBQb3N0cyIsCiAgICAgICAgeWxpbSA9IGMoMCwgbWF4KHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MpICsgNSksCiAgICAgICAgbGFzID0gMiwKICAgICAgICBjb2wgPSAiI0U3NEMzQyIpCgp0ZXh0KHggPSBiYXJwbG90KHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MsIHBsb3QgPSBGQUxTRSksCiAgICAgeSA9IHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MsCiAgICAgbGFiZWxzID0gcm91bmQocG9zdHNfcGVyX3dlZWskYXZnX3Bvc3RzX3Blcl93ZWVrcywgMiksCiAgICAgcG9zID0gMykKYGBgCgojIERlZmluaW5nIHRoZSBpbnRlcnZhbHMgb2YgdGltZSBmb3IgdGhlIGFjYWRlbWljIHllYXIKCmBgYHtyfQppbnRlcnZhbHMgPC0gbGlzdCgKICBpbnRlcnZhbDEgPSBhcy5QT1NJWGN0KGMoIjIwMjItMDgtMzEiLCAiMjAyMi0xMi0xNSIpKSwKICBpbnRlcnZhbDIgPSBhcy5QT1NJWGN0KGMoIjIwMjMtMDEtMDQiLCAiMjAyMy0wNC0wMSIpKSwKICBpbnRlcnZhbDMgPSBhcy5QT1NJWGN0KGMoIjIwMjMtMDQtMTQiLCAiMjAyMy0wNi0xNSIpKQopCmBgYAoKIyBGdW5jdGlvbiB0byBjaGVjayBpZiBhIGRhdGUgZmFsbHMgd2l0aGluIGEgZ2l2ZW4gaW50ZXJ2YWwgb2YgdGltZSBhbmQgYXBwbHkgYXBwcm9wcmlhdGUgQm9vbGVhbgoKYGBge3J9CmNoZWNrX2ludGVydmFsIDwtIGZ1bmN0aW9uKGRhdGUpIHsKICBmb3IgKGkgaW4gMTpsZW5ndGgoaW50ZXJ2YWxzKSkgewogICAgaW50ZXJ2YWxfc3RhcnQgPC0gaW50ZXJ2YWxzW1tpXV1bMV0KICAgIGludGVydmFsX2VuZCA8LSBpbnRlcnZhbHNbW2ldXVsyXQogICAgaWYgKGRhdGUgPj0gaW50ZXJ2YWxfc3RhcnQgJiBkYXRlIDw9IGludGVydmFsX2VuZCkgewogICAgICByZXR1cm4oVFJVRSkKICAgIH0KICB9CiAgcmV0dXJuKEZBTFNFKQp9CmBgYAoKYGBge3J9CmRhdGFfcG9zdHMkYWNhZGVtaWNfeWVhciA8LSBzYXBwbHkoZGF0YV9wb3N0cyRjcmVhdGVkX2F0LCBjaGVja19pbnRlcnZhbCkKcHJpbnQoZGF0YS5mcmFtZShpZCA9IGRhdGFfcG9zdHMkaWQsIGFjYWRlbWljX3llYXIgPSBkYXRhX3Bvc3RzJGFjYWRlbWljX3llYXIpKQpgYGAKCiMgRnVuY3Rpb24gdG8gY291bnQgbnVtYmVyIG9mIHBvc3RzIGFuZCBhdmVyYWdlIHBlciBkYXkgZHVyaW5nIGFjYWRlbWljIHRpbWUgYW5kIHZhY2F0aW9uIHRpbWUKCmBgYHtyfQphbmFseXplX3Bvc3RzIDwtIGZ1bmN0aW9uKGFjYWRlbWljX3llYXJfZmlsdGVyKSB7CiAgIyBGaWx0ZXJpbmcgdGhlIGRhdGEgYmFzZWQgb24gdGhlIGFjYWRlbWljX3llYXJfZmlsdGVyCiAgZmlsdGVyZWRfZGF0YSA8LSBkYXRhX3Bvc3RzICU+JQogICAgZmlsdGVyKGFjYWRlbWljX3llYXIgPT0gYWNhZGVtaWNfeWVhcl9maWx0ZXIpCiAgCiAgIyBDb3VudCBvZiBkYXlzIGZvciBlYWNoIEhFSQogIHVuaXF1ZV9kYXlzIDwtIGZpbHRlcmVkX2RhdGEgJT4lCiAgICBncm91cF9ieShpZCkgJT4lCiAgICBzdW1tYXJpc2UodW5pcXVlX2RheXMgPSBuX2Rpc3RpbmN0KGFzLkRhdGUoY3JlYXRlZF9hdCkpKQogIAogICMgQ291bnQgb2YgcG9zdHMgZm9yIGVhY2ggSEVJCiAgbnVtYmVyX3Bvc3RzX2Jvb2xlYW4gPC0gZmlsdGVyZWRfZGF0YSAlPiUKICAgIGdyb3VwX2J5KGlkKSAlPiUKICAgIHN1bW1hcmlzZShjb3VudCA9IG4oKSkKICAKICAjIE5hbWluZyB0aGUgY29sdW1uIHJlc3BlY3RpbmcgdGhlIHRpbWUgcGVyaW9kCiAgdGltZSA8LSBpZmVsc2UoYWNhZGVtaWNfeWVhcl9maWx0ZXIsICJhY2FkZW1pY190aW1lIiwgInZhY2F0aW9uX3RpbWUiKQogIGNvbHVtbl9uYW1lIDwtIHBhc3RlMCgiYXZnX3Bvc3RzX2luXyIsIHRpbWUpCiAgCiAgIyBDb21iaW5hdGlvbiBvZiBkYXRhIGFuZCBjYWxjdWxhdGlvbiBvZiBhdmVyYWdlIHBvc3RzIHBlciBkYXkKICBjb21iaW5lZF9kYXRhIDwtIGxlZnRfam9pbih1bmlxdWVfZGF5cywgbnVtYmVyX3Bvc3RzX2Jvb2xlYW4sIGJ5ID0gImlkIikKICBjb21iaW5lZF9kYXRhIDwtIGNvbWJpbmVkX2RhdGEgJT4lCiAgICBtdXRhdGUoISFjb2x1bW5fbmFtZSA6PSByb3VuZCgoY291bnQgLyB1bmlxdWVfZGF5cyksIDIpKQogIAogIHByaW50KGNvbWJpbmVkX2RhdGEpCiAgcmV0dXJuKGNvbWJpbmVkX2RhdGEpCn0KYGBgCgpgYGB7cn0KZGF0YV9wb3N0c19hY2FkZW1pYyA8LSBhbmFseXplX3Bvc3RzKFRSVUUpCmRhdGFfcG9zdHNfdmFjYXRpb25zIDwtIGFuYWx5emVfcG9zdHMoRkFMU0UpCmBgYAoKIyBQbG90IGZvciB0aGUgYXZlcmFnZSBudW1iZXIgb2YgcG9zdHMgZHVyaW5nIGFjYWRlbWljIHRpbWUgZm9yIGVhY2ggSEVJCgpgYGB7cn0KYmFycGxvdChkYXRhX3Bvc3RzX2FjYWRlbWljJGF2Z19wb3N0c19pbl9hY2FkZW1pY190aW1lLAogICAgICAgIG5hbWVzLmFyZyA9IGRhdGFfcG9zdHNfYWNhZGVtaWMkaWQsCiAgICAgICAgbWFpbiA9ICJBdmVyYWdlIFBvc3RzIGR1cmluZyBBY2FkZW1pYyBUaW1lIiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJBdmVyYWdlIE51bWJlciBvZiBQb3N0cyIsCiAgICAgICAgeWxpbSA9IGMoMCwgbWF4KGRhdGFfcG9zdHNfYWNhZGVtaWMkYXZnX3Bvc3RzX2luX2FjYWRlbWljX3RpbWUpICsgNSksCiAgICAgICAgbGFzID0gMiwKICAgICAgICBjb2wgPSAiIzM0NDk1RSIpCgp0ZXh0KHggPSBiYXJwbG90KGRhdGFfcG9zdHNfYWNhZGVtaWMkYXZnX3Bvc3RzX2luX2FjYWRlbWljX3RpbWUsIHBsb3QgPSBGQUxTRSksCiAgICAgeSA9IGRhdGFfcG9zdHNfYWNhZGVtaWMkYXZnX3Bvc3RzX2luX2FjYWRlbWljX3RpbWUsCiAgICAgbGFiZWxzID0gcm91bmQoZGF0YV9wb3N0c19hY2FkZW1pYyRhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSwgMiksCiAgICAgcG9zID0gMykKYGBgCgojIFBsb3QgZm9yIHRoZSBhdmVyYWdlIG51bWJlciBvZiBwb3N0cyBkdXJpbmcgdmFjYXRpb24gdGltZSBmb3IgZWFjaCBIRUkKCmBgYHtyfQpiYXJwbG90KGRhdGFfcG9zdHNfdmFjYXRpb25zJGF2Z19wb3N0c19pbl92YWNhdGlvbl90aW1lLAogICAgICAgIG5hbWVzLmFyZyA9IGRhdGFfcG9zdHNfdmFjYXRpb25zJGlkLAogICAgICAgIG1haW4gPSAiQXZlcmFnZSBQb3N0cyBkdXJpbmcgVmFjYXRpb24gVGltZSIsCiAgICAgICAgeGxhYiA9ICJIRUkiLAogICAgICAgIHlsYWIgPSAiQXZlcmFnZSBOdW1iZXIgb2YgUG9zdHMiLAogICAgICAgIHlsaW0gPSBjKDAsIG1heChkYXRhX3Bvc3RzX3ZhY2F0aW9ucyRhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSkgKyA1KSwKICAgICAgICBsYXMgPSAyLAogICAgICAgIGNvbCA9ICIjRDM1NDAwIikKCnRleHQoeCA9IGJhcnBsb3QoZGF0YV9wb3N0c192YWNhdGlvbnMkYXZnX3Bvc3RzX2luX3ZhY2F0aW9uX3RpbWUsIHBsb3QgPSBGQUxTRSksCiAgICAgeSA9IGRhdGFfcG9zdHNfdmFjYXRpb25zJGF2Z19wb3N0c19pbl92YWNhdGlvbl90aW1lLAogICAgIGxhYmVscyA9IHJvdW5kKGRhdGFfcG9zdHNfdmFjYXRpb25zJGF2Z19wb3N0c19pbl92YWNhdGlvbl90aW1lLCAyKSwKICAgICBwb3MgPSAzKQpgYGAKCiMgRGF0YSBwcmVwYXJhdGlvbiBmb3IgZGF0ZXMgCgpgYGB7cn0KIyBDcmVhdGluZyBuZXcgdGFibGUgdGhhdCBjb250YWlucyBhIG5ldyBjb2x1bW4gZm9yIHRoZSBkYXkgb2YgdGhlIHdlZWsKZGF0YV9wb3N0c19kYXlzIDwtIGRhdGFfcG9zdHMgJT4lCiAgbXV0YXRlKGRheV9vZl93ZWVrID0gd2Vla2RheXMoY3JlYXRlZF9hdCkpCgojIFNlbGVjdGluZyBvbmx5IHRoZSBpZCwgY3JlYXRlZF9hdCwgYW5kIGRheV9vZl93ZWVrIGNvbHVtbnMgZm9yIHRoZSBuZXcgdGFibGUKZGF0YV9wb3N0c19kYXlzIDwtIGRhdGFfcG9zdHNfZGF5cyAlPiUKICBzZWxlY3QoaWQsIGNyZWF0ZWRfYXQsIGRheV9vZl93ZWVrKQoKIyBDcmVhdGUgY29sdW1uIGhvdXIgZnJvbSBjcmVhdGVkX2F0CmRhdGFfcG9zdHNfZGF5cyRjcmVhdGVkX2hvdXIgPC0gYXMubnVtZXJpYyhmb3JtYXQoZGF0YV9wb3N0c19kYXlzJGNyZWF0ZWRfYXQsICIlSCIpKQoKcHJpbnQoZGF0YV9wb3N0c19kYXlzKQpgYGAKCmBgYHtyfQojIEdyb3VwaW5nIGJ5IGlkIGFuZCBkYXlfb2Zfd2VlaywgdGhlbiBjb3VudGluZyB0aGUgbnVtYmVyIG9mIHBvc3RzCm51bWJlcl9wb3N0c19kYXlzIDwtIGRhdGFfcG9zdHNfZGF5cyAlPiUKICBncm91cF9ieShpZCwgZGF5X29mX3dlZWspICU+JQogIHN1bW1hcmlzZShjb3VudCA9IG4oKSkKCiMgR3JvdXBpbmcgYnkgaWQsIGRheV9vZl93ZWVrIGFuZCBkYXkgY3JlYXRlZCBhdCwgdGhlbiBjb3VudGluZyB0aCBlbnVtYmVyIG9mIHR3ZWV0cwpudW1iZXJfcG9zdHNfcGVyX2RheSA8LSBkYXRhX3Bvc3RzX2RheXMgJT4lCiAgICBtdXRhdGUoY3JlYXRlZF9kYXRlID0gYXMuRGF0ZShjcmVhdGVkX2F0KSkgJT4lIAogICAgZ3JvdXBfYnkoaWQsIGRheV9vZl93ZWVrLCBjcmVhdGVkX2RhdGUpICU+JQogICAgc3VtbWFyaXplKGNvdW50ID0gbigpKQoKIyBGaW5kaW5nIGZvciBlYWNoIEhFSSB0aGUgYXZlcmFnZSBjb3VudCBvZiBwb3N0cyBwZXIgZGF5CmF2ZXJhZ2VfbnVtYmVyX3Bvc3RzX3Blcl9kYXkgPC0gbnVtYmVyX3Bvc3RzX3Blcl9kYXkgJT4lCiAgZ3JvdXBfYnkoaWQsIGRheV9vZl93ZWVrKSAlPiUKICBzdW1tYXJpc2UoYXZlcmFnZV9jb3VudCA9IHJvdW5kKG1lYW4oY291bnQpLCAyKSkKCnByaW50KG51bWJlcl9wb3N0c19kYXlzKQpgYGAKCiMgSGlnaGVzdCBhbmQgbG93ZXN0IHBvc3RzCgpgYGB7cn0KIyBGaW5kaW5nIHRoZSBIRUkgd2l0aCB0aGUgbG93ZXN0IGNvdW50IG9mIHBvc3RzIHBlciBkYXkKbG93ZXN0X2NvdW50IDwtIG51bWJlcl9wb3N0c19kYXlzICU+JQogIGdyb3VwX2J5KGRheV9vZl93ZWVrKSAlPiUKICBzbGljZV9taW4ob3JkZXJfYnkgPSBjb3VudCkgJT4lCiAgc2VsZWN0KGRheV9vZl93ZWVrLCBpZCwgY291bnQpCgojIEZpbmRpbmcgdGhlIEhFSSB3aXRoIHRoZSBoaWdoZXN0IGNvdW50IG9mIHBvc3RzIHBlciBkYXkKaGlnaGVzdF9jb3VudCA8LSBudW1iZXJfcG9zdHNfZGF5cyAlPiUKICBncm91cF9ieShkYXlfb2Zfd2VlaykgJT4lCiAgc2xpY2VfbWF4KG9yZGVyX2J5ID0gY291bnQpICU+JQogIHNlbGVjdChkYXlfb2Zfd2VlaywgaWQsIGNvdW50KQoKIyBDb21iaW5lIHRoZSByZXN1bHRzCmhpZ2hfbG93X0hFSSA8LSBiaW5kX3Jvd3MobG93ZXN0X2NvdW50LCBoaWdoZXN0X2NvdW50KSAlPiUKICBhcnJhbmdlKGRheV9vZl93ZWVrKQoKcHJpbnQoaGlnaF9sb3dfSEVJKQpgYGAKCiMgUGxvdCBmb3IgdGhlIGhpZ2hlc3QgYW5kIGxvd2VzdCBjb3VudCBvZiBwb3N0cyBwZXIgZGF5IGZvciBlYWNoIGRheSBvZiB0aGUgd2VlawoKYGBge3J9CmdncGxvdChoaWdoX2xvd19IRUksIGFlcyh4ID0gZGF5X29mX3dlZWssIHkgPSBjb3VudCwgZmlsbCA9IGlkKSkgKwogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsKICBnZW9tX3RleHQoYWVzKGxhYmVsID0gY291bnQpLAogICAgICAgICAgICBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC45KSwKICAgICAgICAgICAgdmp1c3QgPSAtMC41LAogICAgICAgICAgICBzaXplID0gMykgKwogIGxhYnModGl0bGUgPSAiSGlnaGVzdCBhbmQgTG93ZXN0IENvdW50IG9mIFBvc3RzIHBlciBEYXkgZm9yIEVhY2ggRGF5IG9mIHRoZSBXZWVrIiwKICAgICAgIHggPSAiRGF5IG9mIHRoZSBXZWVrIiwgeSA9ICJDb3VudCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSByYWluYm93KGxlbmd0aCh1bmlxdWUoaGlnaF9sb3dfSEVJJGlkKSkpKSArCiAgdGhlbWVfbWluaW1hbCgpICsKICB0aGVtZShsZWdlbmQudGl0bGUgPSBlbGVtZW50X2JsYW5rKCkpCmBgYAoKIyBBdmVyYWdlIG9mIHBvc3RzCgpgYGB7cn0KIyBGaW5kaW5nIHRoZSBIRUkgd2l0aCBsb3dlc3QgYW5kIGhpZ2hlc3QgYXZlcmFnZWQgY291bnQgb2YgcG9zdHMgcGVyIGRheQpoaWdoX2xvd19hdmVyYWdlX0hFSXMgPC0gYXZlcmFnZV9udW1iZXJfcG9zdHNfcGVyX2RheSAlPiUKICBncm91cF9ieShkYXlfb2Zfd2VlaykgJT4lCiAgZmlsdGVyKGF2ZXJhZ2VfY291bnQgPT0gbWF4KGF2ZXJhZ2VfY291bnQpIHwgYXZlcmFnZV9jb3VudCA9PSBtaW4oYXZlcmFnZV9jb3VudCkpICU+JQogIGFycmFuZ2UoZGF5X29mX3dlZWssIGlmZWxzZShhdmVyYWdlX2NvdW50ID09IG1pbihhdmVyYWdlX2NvdW50KSwgYXZlcmFnZV9jb3VudCwgLWF2ZXJhZ2VfY291bnQpKQoKcHJpbnQoaGlnaF9sb3dfYXZlcmFnZV9IRUlzKQpgYGAKCiMgUGxvdCBmb3IgdGhlIGhpZ2hlc3QgYW5kIGxvd2VzdCBhdmVyYWdlIGNvdW50IG9mIHBvc3RzIHBlciBkYXkgZm9yIGVhY2ggZGF5IG9mIHRoZSB3ZWVrCgpgYGB7cn0KZ2dwbG90KGhpZ2hfbG93X2F2ZXJhZ2VfSEVJcywgYWVzKHggPSBkYXlfb2Zfd2VlaywgeSA9IGF2ZXJhZ2VfY291bnQsIGZpbGwgPSBpZCkpICsKICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5IiwgcG9zaXRpb24gPSAiZG9kZ2UiKSArCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbCA9IHJvdW5kKGF2ZXJhZ2VfY291bnQsIDIpKSwKICAgICAgICAgICAgcG9zaXRpb24gPSBwb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNyksCiAgICAgICAgICAgIHZqdXN0ID0gLTAuNSwKICAgICAgICAgICAgc2l6ZSA9IDMpICsKICBsYWJzKHRpdGxlID0gIkhpZ2hlc3QgYW5kIExvd2VzdCBBdmVyYWdlIENvdW50IG9mIFBvc3RzIHBlciBEYXkgZm9yIEVhY2ggRGF5IG9mIHRoZSBXZWVrIiwKICAgICAgIHggPSAiRGF5IG9mIHRoZSBXZWVrIiwgeSA9ICJBdmVyYWdlIENvdW50IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHJhaW5ib3cobGVuZ3RoKHVuaXF1ZShoaWdoX2xvd19IRUkkaWQpKSkpICsKICB0aGVtZV9taW5pbWFsKCkgKwogIHRoZW1lKGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfYmxhbmsoKSkKYGBgCgojIEZhdm91cml0ZSBob3VyIGFuZCBkYXkKCmBgYHtyfQpmYXZvdXJpdGVfZGF5X2hlaSA8LSBudW1iZXJfcG9zdHNfZGF5cyAlPiUKICBncm91cF9ieShpZCkgJT4lCiAgdG9wX24oMSwgY291bnQpICU+JQogIGFycmFuZ2UoaWQpCgpwcmludChmYXZvdXJpdGVfZGF5X2hlaSkKYGBgCgpgYGB7cn0KIyBDYWxjdWxhdGluZyB0aGUgbnVtYmVyIG9mIHBvc3RzIHBlciBob3VyCm51bWJlcl9wb3N0c19ob3VycyA8LSBkYXRhX3Bvc3RzX2RheXMgJT4lCiAgZ3JvdXBfYnkoaWQsIGNyZWF0ZWRfaG91cikgJT4lCiAgc3VtbWFyaXNlKGNvdW50ID0gbigpKSAlPiUKICB1bmdyb3VwKCkKCiMgSWRlbnRpZnlpbmcgdGhlIGZhdm9yaXRlIGhvdXIgZm9yIGVhY2ggSEVJCmZhdm91cml0ZV9ob3VyX2hlaSA8LSBudW1iZXJfcG9zdHNfaG91cnMgJT4lCiAgZ3JvdXBfYnkoaWQpICU+JQogIHRvcF9uKDEsIGNvdW50KSAlPiUKICBhcnJhbmdlKGlkKQoKIyBBZGRpbmcgbmV3IGNvbHVtbnMgdG8gY2F0ZWdvcml6ZSB0aGUgaG91cnMgYW5kIGFzc2lnbiBudW1lcmljYWwgdmFsdWUgZm9yIGNsdXN0ZXJpbmcKZmF2b3VyaXRlX2hvdXJfaGVpIDwtIGZhdm91cml0ZV9ob3VyX2hlaSAlPiUKICBtdXRhdGUodGltZV9vZl9kYXkgPSBjYXNlX3doZW4oCiAgICBjcmVhdGVkX2hvdXIgPj0gNSAmIGNyZWF0ZWRfaG91ciA8IDEyIH4gIk1vcm5pbmciLAogICAgY3JlYXRlZF9ob3VyID49IDEyICYgY3JlYXRlZF9ob3VyIDwgMTggfiAiQWZ0ZXJub29uIiwKICAgIFRSVUUgfiAiTmlnaHQiCiAgKSwKICB0aW1lX29mX2RheV92YWx1ZSA9IGNhc2Vfd2hlbigKICAgIGNyZWF0ZWRfaG91ciA+PSA1ICYgY3JlYXRlZF9ob3VyIDwgMTIgfiAxLAogICAgY3JlYXRlZF9ob3VyID49IDEyICYgY3JlYXRlZF9ob3VyIDwgMTggfiAzLAogICAgVFJVRSB+IDUKICApKQoKcHJpbnQoZmF2b3VyaXRlX2hvdXJfaGVpKQpgYGAKCiMgSGVhdG1hcHMKCiMgRnVuY3Rpb24gdG8gcGxvdCBoZWF0bWFwIGZvciB2YXJpb3VzIEhFSXMKCmBgYHtyfQpoZWF0bWFwX21ha2VyIDwtIGZ1bmN0aW9uKHRhcmdldF9pZCl7CiAgIyBGaWx0ZXJpbmcgZGF0YSBmb3IgdGhlIHNwZWNpZmljIEhFSQogIHRhcmdldF9kYXRhIDwtIGRhdGFfcG9zdHNfZGF5cyAlPiUKICAgIGZpbHRlcihpZCA9PSB0YXJnZXRfaWQpCiAgCiAgIyBHcm91cGluZyBieSBkYXkgb2YgdGhlIHdlZWsgYW5kIGhvdXIsIGFuZCBjb3VudGluZyB0aGUgbnVtYmVyIG9mIHR3ZWV0cwogIHR3ZWV0X2NvdW50cyA8LSB0YXJnZXRfZGF0YSAlPiUKICAgIGdyb3VwX2J5KGRheV9vZl93ZWVrLCBjcmVhdGVkX2hvdXIpICU+JQogICAgc3VtbWFyaXNlKG51bV9wb3N0cyA9IG4oKSkKICAKICAjIFBsb3R0aW5nIGhlYXRtYXAKICBnZ3Bsb3QodHdlZXRfY291bnRzLCBhZXMoeCA9IGRheV9vZl93ZWVrLCB5ID0gY3JlYXRlZF9ob3VyLCBmaWxsID0gbnVtX3Bvc3RzKSkgKwogICAgZ2VvbV90aWxlKCkgKwogICAgc2NhbGVfZmlsbF9ncmFkaWVudChsb3cgPSAid2hpdGUiLCBoaWdoID0gImJsdWUiKSArCiAgICBsYWJzKHRpdGxlID0gcGFzdGUoIlBvc3QgSGVhdG1hcCBmb3IiLCB0YXJnZXRfaWQpLAogICAgICAgICB4ID0gIkRheSBvZiB0aGUgd2VlayIsCiAgICAgICAgIHkgPSAiSG91ciBvZiB0aGUgZGF5IikKfQpgYGAKCiMgUGxvdCBvZiBoZWF0bWFwIGZvciBlYWNoIEhFSQoKYGBge3J9CmhlYXRtYXBfbWFrZXIoImR1a2UuY3N2IikKaGVhdG1hcF9tYWtlcigiZXBmbC5jc3YiKQpoZWF0bWFwX21ha2VyKCJnb2UuY3N2IikKaGVhdG1hcF9tYWtlcigiaGFydmFyZC5jc3YiKQpoZWF0bWFwX21ha2VyKCJsZWljZXN0ZXIuY3N2IikKaGVhdG1hcF9tYWtlcigibWFuY2hlc3Rlci5jc3YiKQpoZWF0bWFwX21ha2VyKCJtaXQuY3N2IikKaGVhdG1hcF9tYWtlcigic2IuY3N2IikKaGVhdG1hcF9tYWtlcigic3RhbmZvcmQuY3N2IikKaGVhdG1hcF9tYWtlcigidHJpbml0eS5jc3YiKQpoZWF0bWFwX21ha2VyKCJ3di5jc3YiKQpoZWF0bWFwX21ha2VyKCJ5YWxlLmNzdiIpCmBgYAoKIyBIYXNodGFncwoKYGBge3J9CiMgVHJhbnNmb3JtaW5nIGVtcHR5IHN0cmluZ3MgaW50byBOQQpkYXRhX3Bvc3RzJGhhc2h0YWdzW2RhdGFfcG9zdHMkaGFzaHRhZ3MgPT0gIiJdIDwtIE5BCgojIFRhYmxlIHdpdGggbnVtYmVyIG9mIHVuaXF1ZSBoYXNodGFncyBhbmQgcGVyY2VudGFnZSBvZiB1c2FnZQpoYXNodGFncyA8LSBkYXRhX3Bvc3RzICU+JQogICAgICAgICAgICAgICAgZ3JvdXBfYnkoaWQpICU+JQogICAgICAgICAgICAgICAgc3VtbWFyaXNlKGNvdW50ID0gbigpLAogICAgICAgICAgICAgICAgICAgICAgICAgIG5hID0gc3VtKGlzLm5hKGhhc2h0YWdzKSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgdW5pcXVlX2hhc2h0YWdzID0gbGVuZ3RoKHVuaXF1ZShoYXNodGFncykpLAogICAgICAgICAgICAgICAgICAgICAgICAgIGhhc2h0YWdfcGVyY2VudGFnZSA9IHJvdW5kKCgoY291bnQgLSBuYSkgLyBjb3VudCAqIDEwMCksIDIpKQoKcHJpbnQoaGFzaHRhZ3MpCmBgYAoKIyBQbG90IGZvciB0aGUgY291bnQgb2YgdW5pcXVlIGhhc2h0YWdzIGZvciBlYWNoIEhFSQoKYGBge3J9CmJhcnBsb3QoaGFzaHRhZ3MkdW5pcXVlX2hhc2h0YWdzLAogICAgICAgIG5hbWVzLmFyZyA9IGhhc2h0YWdzJGlkLAogICAgICAgIG1haW4gPSAiVW5pcXVlIEhhc2h0YWdzIGZvciBFYWNoIEhFSSIsCiAgICAgICAgeGxhYiA9ICJIRUkiLAogICAgICAgIHlsYWIgPSAiQ291bnQgb2YgVW5pcXVlIEhhc2h0YWdzIiwKICAgICAgICB5bGltID0gYygwLCBtYXgoaGFzaHRhZ3MkdW5pcXVlX2hhc2h0YWdzKSArIDUwKSwKICAgICAgICBsYXMgPSAyLAogICAgICAgIGNvbD0gIiMxNkEwODUiKQoKdGV4dCh4ID0gYmFycGxvdChoYXNodGFncyR1bmlxdWVfaGFzaHRhZ3MsIHBsb3QgPSBGQUxTRSksCiAgICAgeSA9IGhhc2h0YWdzJHVuaXF1ZV9oYXNodGFncywKICAgICBsYWJlbHMgPSByb3VuZChoYXNodGFncyR1bmlxdWVfaGFzaHRhZ3MsIDIpLAogICAgIHBvcyA9IDMpCmBgYAoKIyBQbG90IGZvciB0aGUgdXNhZ2Ugb2YgaGFzaHRhZyBmb3IgZWFjaCBIRUkKCmBgYHtyfQpiYXJwbG90KGhhc2h0YWdzJGhhc2h0YWdfcGVyY2VudGFnZSwKICAgICAgICBuYW1lcy5hcmcgPSBoYXNodGFncyRpZCwKICAgICAgICBtYWluID0gIkhhc2h0YWdzIFBlcmNlbnRhZ2UgZm9yIEVhY2ggSEVJIiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJIYXNodGFncyBQZXJjZW50YWdlIiwKICAgICAgICB5bGltID0gYygwLCBtYXgoaGFzaHRhZ3MkaGFzaHRhZ19wZXJjZW50YWdlKSArIDMwKSwKICAgICAgICBsYXMgPSAyLAogICAgICAgIGNvbD0gIiNGMUM0MEYiKQoKdGV4dCh4ID0gYmFycGxvdChoYXNodGFncyRoYXNodGFnX3BlcmNlbnRhZ2UsIHBsb3QgPSBGQUxTRSksCiAgICAgeSA9IGhhc2h0YWdzJGhhc2h0YWdfcGVyY2VudGFnZSwKICAgICBsYWJlbHMgPSByb3VuZChoYXNodGFncyRoYXNodGFnX3BlcmNlbnRhZ2UsIDIpLAogICAgIHBvcyA9IDMpCmBgYAoKIyBVUkwgdXNhZ2UKCmBgYHtyfQojIFRyYW5zZm9ybWluZyBlbXB0eSBzdHJpbmdzIGludG8gTkEKZGF0YV9wb3N0cyR1cmxzW2RhdGFfcG9zdHMkdXJscyA9PSAiIl0gPC0gTkEKCiMgVGFibGUgd2l0aCBudW1iZXIgb2YgcG9zdCwgbnVtYmVyIG9mIE5BIGFuZCB1cmwgcGVyY2VudGFnZSBvZiB1c2FnZQp1cmxfdXNhZ2UgPC0gZGF0YV9wb3N0cyAlPiUKICAgICAgICAgICAgICAgIGdyb3VwX2J5KGlkKSAlPiUKICAgICAgICAgICAgICAgIHN1bW1hcmlzZShjb3VudCA9IG4oKSwKICAgICAgICAgICAgICAgICAgICAgICAgICBuYSA9IHN1bShpcy5uYSh1cmxzKSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgdXJsX3BlcmNlbnRhZ2UgPSByb3VuZCgoKGNvdW50IC0gbmEpIC8gY291bnQgKiAxMDApLCAyKSkKCnByaW50KHVybF91c2FnZSkKYGBgCgojIFBsb3QgZm9yIHRoZSB1c2FnZSBvZiBoYXNodGFnIGZvciBlYWNoIEhFSQoKYGBge3J9CmJhcnBsb3QodXJsX3VzYWdlJHVybF9wZXJjZW50YWdlLAogICAgICAgIG5hbWVzLmFyZyA9IHVybF91c2FnZSRpZCwKICAgICAgICBtYWluID0gIlVybHMgUGVyY2VudGFnZSBmb3IgRWFjaCBIRUkiLAogICAgICAgIHhsYWIgPSAiSEVJIiwKICAgICAgICB5bGFiID0gIlVybHMgUGVyY2VudGFnZSIsCiAgICAgICAgeWxpbSA9IGMoMCwgbWF4KHVybF91c2FnZSR1cmxfcGVyY2VudGFnZSkgKyAxMCksCiAgICAgICAgbGFzID0gMiwKICAgICAgICBjb2w9ICIjOEU0NEFEIikKCnRleHQoeCA9IGJhcnBsb3QodXJsX3VzYWdlJHVybF9wZXJjZW50YWdlLCBwbG90ID0gRkFMU0UpLAogICAgIHkgPSB1cmxfdXNhZ2UkdXJsX3BlcmNlbnRhZ2UsCiAgICAgbGFiZWxzID0gcm91bmQodXJsX3VzYWdlJHVybF9wZXJjZW50YWdlLCAyKSwKICAgICBwb3MgPSAzKQpgYGAKCiMgVGV4dAoKYGBge3J9CmRhdGFfcG9zdHNfY29udGVudCA8LSBkYXRhX3Bvc3RzICU+JQogICAgICAgICAgICBzZWxlY3QoaWQsIHRleHQpCgojIENvdW50aW5nIG51bWJlciBvZiB3b3JkcwpkYXRhX3Bvc3RzX2NvbnRlbnQgPC0gZGF0YV9wb3N0c19jb250ZW50ICU+JQogIG11dGF0ZShudW1fd29yZHMgPSBsZW5ndGhzKHN0cnNwbGl0KHRleHQsICJcXHMrIikpKQoKIyBHcm91cGluZyBieSBIRUkgYW5kIGNhbGN1bGF0ZSBhdmVyYWdlLCBtaW5pbXVtLCBhbmQgbWF4aW11bSB2YWx1ZXMgb2YgbnVtYmVyIG9mIHdvcmRzCmRhdGFfcG9zdHNfY29udGVudF9tZXRyaWNzIDwtIGRhdGFfcG9zdHNfY29udGVudCAlPiUKICBncm91cF9ieShpZCkgJT4lCiAgc3VtbWFyaXNlKGF2ZXJhZ2VfbnVtX3dvcmRzID0gbWVhbihudW1fd29yZHMpLAogICAgICAgICAgICBtaW5fbnVtX3dvcmRzID0gbWluKG51bV93b3JkcyksCiAgICAgICAgICAgIG1heF9udW1fd29yZHMgPSBtYXgobnVtX3dvcmRzKSkKcHJpbnQoZGF0YV9wb3N0c19jb250ZW50X21ldHJpY3MpCmBgYAoKIyBQbG90IGZvciB0aGUgYXZlcmFnZSwgbWF4aW11bSBhbmQgbWluaW11bSB2YWx1ZXMgb2Ygd29yZHMgZm9yIGVhY2ggSEVJCgpgYGB7cn0KZ2dwbG90KGRhdGFfcG9zdHNfY29udGVudF9tZXRyaWNzLCBhZXMoeCA9IGlkLCB5ID0gYXZlcmFnZV9udW1fd29yZHMpKSArCiAgZ2VvbV9wb2ludChhZXMoY29sb3IgPSAiQXZlcmFnZSIpKSArCiAgZ2VvbV9lcnJvcmJhcihhZXMoeW1pbiA9IG1pbl9udW1fd29yZHMsIHltYXggPSBtYXhfbnVtX3dvcmRzLCBjb2xvciA9ICJSYW5nZSIpLCB3aWR0aCA9IDAuMikgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJBdmVyYWdlIiA9ICIjMTk3NkQyIiwgIlJhbmdlIiA9ICIjRUY1MzUwIikpICsKICBsYWJzKHRpdGxlID0gIldvcmQgQ291bnQgU3VtbWFyeSBieSBIRUkiLAogICAgICAgeCA9ICJIRUkiLAogICAgICAgeSA9ICJOdW1iZXIgb2YgV29yZHMiLAogICAgICAgY29sb3IgPSAiTWV0cmljIikgKwogIHRoZW1lX21pbmltYWwoKSArCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0ID0gMSkpCmBgYAoKIyBQb3N0cyBDbGFzc2lmaWNhdGlvbgoKYGBge3J9CmNsZWFudXAgPC0gZnVuY3Rpb24oZG9jcywgc3BlYy53b3Jkcz1OVUxMKXsKICAjIGxvd2VyY2FzZQogIGRvY3MgPC0gdG1fbWFwKGRvY3MsIGNvbnRlbnRfdHJhbnNmb3JtZXIodG9sb3dlcikpCiAgIyBybSBudW1iZXJzCiAgZG9jcyA8LSB0bV9tYXAoZG9jcywgcmVtb3ZlTnVtYmVycykKICAjIHJtIGVuZ2xpc2ggY29tbW9uIHN0b3BXb3JkcwogIGRvY3MgPC0gdG1fbWFwKGRvY3MsIHJlbW92ZVdvcmRzLCBzdG9wd29yZHMoImVuZ2xpc2giKSkKICAjIGlmIHN0b3B3b3JkcyBhcmUgc3BlY2lmaWVkIGFzIGEgY2hhcmFjdGVyIHZlY3RvcgogIGlmKCFpcy5udWxsKHNwZWMud29yZHMpKQogICAgZG9jcyA8LSB0bV9tYXAoZG9jcywgcmVtb3ZlV29yZHMsIHNwZWMud29yZHMpCiAgIyBybSBwdW5jdHVhdGlvbnMKICBkb2NzIDwtIHRtX21hcChkb2NzLCByZW1vdmVQdW5jdHVhdGlvbikKICAjIHJtIGV4dHJhIHdoaXRlIHNwYWNlcwogIGRvY3MgPC0gdG1fbWFwKGRvY3MsIHN0cmlwV2hpdGVzcGFjZSkKICAjIGxlbW1hdGl6aW5nIHRleHQKICBkb2NzIDwtIHRtX21hcChkb2NzLCBsZW1tYXRpemVfd29yZHMpCiAgCiAgZG9jcwp9CmBgYAoKYGBge3J9CnBhaXJpbmcgPC0gZnVuY3Rpb24od29yZF9kaWN0LCBwYWlyZWRfd29yZHMpIHsKICBpZiAobnJvdyh3b3JkX2RpY3QpICE9IGxlbmd0aChwYWlyZWRfd29yZHMpKSB7CiAgICBzdG9wKCJUaGUgbnVtYmVyIG9mIHJvd3MgaW4gd29yZF9kaWN0IGFuZCB0aGUgbGVuZ3RoIG9mIHBhaXJlZF93b3JkcyBzaG91bGQgYmUgZXF1YWwuIikKICB9CiAgCiAgd29yZF9kaWN0IDwtIHRpYmJsZTo6dGliYmxlKHdvcmQgPSB3b3JkX2RpY3QkdmFsdWUpCiAgcmVzdWx0IDwtIHRpYmJsZTo6dGliYmxlKHdvcmQgPSB3b3JkX2RpY3Qkd29yZCwgY2F0ZWdvcnkgPSBwYWlyZWRfd29yZHMpCiAgcmV0dXJuKHJlc3VsdCkKfQpgYGAKCmBgYHtyfQpjbGFzc2lmeV90ZXh0IDwtIGZ1bmN0aW9uKHRleHQsIHdvcmRfcGFpcikgewogICMgVG9rZW5pemluZyB0aGUgdGV4dCBhbmQgY29udmVydGluZyB0byBsb3dlcmNhc2UKICB3b3JkcyA8LSB0b2xvd2VyKHVubGlzdChzdHJzcGxpdCh0ZXh0LCAiXFxXKyIpKSkKICAKICAjIEZpbmRpbmcgdGhlIGZyZXF1ZW50IHdvcmRzIGluIHRoZSB0ZXh0CiAgZnJlcV93b3JkcyA8LSB3b3Jkc1t3b3JkcyAlaW4lIHdvcmRfcGFpciR3b3JkXQogIAogIGlmIChsZW5ndGgoZnJlcV93b3JkcykgPT0gMCkgewogICAgcmV0dXJuKCJVbmtub3duIikgICMgUmV0dXJuaW5nIFVua25vd24gaWYgbm8gZnJlcXVlbnQgd29yZHMgYXJlIGZvdW5kCiAgfQogIAogICMgR2V0dGluZyB0aGUgY29ycmVzcG9uZGluZyBjYXRlZ29yaWVzCiAgY2F0ZWdvcmllcyA8LSB3b3JkX3BhaXIkY2F0ZWdvcnlbd29yZF9wYWlyJHdvcmQgJWluJSBmcmVxX3dvcmRzXQogIAogICMgU29ydGluZyBjYXRlZ29yaWVzIGJ5IGZyZXF1ZW5jeSBpbiBkZXNjZW5kaW5nIG9yZGVyIGFuZCByZXR1cm4gdGhlIG1vc3QgZnJlcXVlbnQgb25lCiAgcmV0dXJuKG5hbWVzKHNvcnQodGFibGUoY2F0ZWdvcmllcyksIGRlY3JlYXNpbmcgPSBUUlVFKSlbMV0pCn0KYGBgCgpgYGB7cn0KY29ycHVzX21ha2VyIDwtIGZ1bmN0aW9uKHRleHQpIHsKICB0ZXh0cyA8LSB0ZXh0JHRleHQKICB2YyA8LSBWZWN0b3JTb3VyY2UodGV4dHMpCiAgY29ycHVzIDwtIENvcnB1cyh2YykKICAKICByZXR1cm4oY29ycHVzKQp9CmBgYAoKYGBge3J9CmZyZXFfdGVybXMgPC0gZnVuY3Rpb24oY2xlYW5fdGV4dCwgbnVtYmVyKSB7CiAgZHRtIDwtIERvY3VtZW50VGVybU1hdHJpeChjbGVhbl90ZXh0KQogIGR0bS50ZmlkZiA8LSB3ZWlnaHRUZklkZihkdG0pCiAgCiAgbWRmIDwtIGFzX3RpYmJsZShhcy5tYXRyaXgoZHRtLnRmaWRmKSkKICAKICBtZGYuZnJlcSA8LSBtZGYgJT4lCiAgc2VsZWN0KGZpbmRGcmVxVGVybXMoZHRtLCBudW1iZXIpKSAlPiUKICBzdW1tYXJpc2VfYWxsKHN1bSkgJT4lCiAgZ2F0aGVyKCkgJT4lCiAgYXJyYW5nZShkZXNjKHZhbHVlKSkKCiAgbWRmLmZyZXEka2V5IDwtIAogICAgZmFjdG9yKG1kZi5mcmVxJGtleSwKICAgICAgICAgICBsZXZlbHMgPSBtZGYuZnJlcSRrZXlbb3JkZXIobWRmLmZyZXEkdmFsdWUpXSkKICAKICB3b3JkX2RpY3Rpb25hcnkgPC0gYXNfdGliYmxlKG1kZi5mcmVxJGtleSkKICAKICBnZ3Bsb3QobWRmLmZyZXEsIGFlcyh4PWtleSwgeT12YWx1ZSkpICsKICAgIGdlb21fYmFyKHN0YXQ9ImlkZW50aXR5IikgKwogICAgbGFicyh4PSJ0ZXJtcyIsIHk9ImZyZXEiKSArIGNvb3JkX2ZsaXAoKQogIAogIHJldHVybih3b3JkX2RpY3Rpb25hcnkpCn0KYGBgCgojIER1a2UKCmBgYHtyfQpkdWtlX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gImR1a2UuY3N2IikKZHVrZV9jb3JwdXMgPC0gY29ycHVzX21ha2VyKGR1a2VfdGV4dCkKCmR1a2VfY2xlYW5fdGV4dCA8LWNsZWFudXAoZHVrZV9jb3JwdXMsIGMoIm5ldyIsICJ3aWxsIiwgImNoYW5nZSIsICJub3J0aCIsICJjYW4iLCAiZmlyc3QiLCAieWVhciIsICJjYXJvbGluYSIsICJ5ZWFycyIsICJzdW1tZXIiLCAibWVldCIsICJvbmUiKSkKCmR1a2Vfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXMoZHVrZV9jbGVhbl90ZXh0LCAzMCkKCmR1a2Vfd29yZF9jYXRlZ29yeSA8LSBjKCJJbWFnZSIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIkltYWdlIiwgIkVuZ2FnZW1lbnQiLCAiUmVzZWFyY2giLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIsICJTb2NpZXR5IiwgIkVkdWNhdGlvbiIsICJFbmdhZ2VtZW50IikKCmR1a2Vfd29yZF9wYWlyIDwtIHBhaXJpbmcoZHVrZV93b3JkX2RpY3Rpb25hcnksIGR1a2Vfd29yZF9jYXRlZ29yeSkKCnByaW50KGR1a2Vfd29yZF9wYWlyKQoKIyBBcHBseWluZyB0aGUgY2xhc3NpZnlfdGV4dCBmdW5jdGlvbiB0byBlYWNoIHRleHQKZHVrZV90ZXh0IDwtIGR1a2VfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gZHVrZV93b3JkX3BhaXIpKQoKcHJpbnQoZHVrZV90ZXh0KQpgYGAKCiMgZXBmbAoKYGBge3J9CmVwZmxfdGV4dCA8LSBzdWJzZXQoZGF0YV9wb3N0c19jb250ZW50LCBpZCA9PSAiZXBmbC5jc3YiKQplcGZsX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoZXBmbF90ZXh0KQoKZXBmbF9jbGVhbl90ZXh0IDwtIGNsZWFudXAoZXBmbF9jb3JwdXMsIGMoIm5ldyIsICJhbXAiLCAiY2FuIiwgIndpbGwiLCAib25lIiwgIm5vdyIsICLigJMiLCAidmlhIiwgInBvcnRyYWl0IikpCgplcGZsX3dvcmRfZGljdGlvbmFyeSA8LSBmcmVxX3Rlcm1zKGVwZmxfY2xlYW5fdGV4dCwgMzApCgplcGZsX3dvcmRfY2F0ZWdvcnkgPC0gYygiSW1hZ2UiLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiUmVzZWFyY2giLCAiUmVzZWFyY2giLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIsIE5BLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIpCgplcGZsX3dvcmRfcGFpciA8LSBwYWlyaW5nKGVwZmxfd29yZF9kaWN0aW9uYXJ5LCBlcGZsX3dvcmRfY2F0ZWdvcnkpCgpwcmludChlcGZsX3dvcmRfcGFpcikKCiMgQXBwbHlpbmcgdGhlIGNsYXNzaWZ5X3RleHQgZnVuY3Rpb24gdG8gZWFjaCB0ZXh0CmVwZmxfdGV4dCA8LSBlcGZsX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IGVwZmxfd29yZF9wYWlyKSkKCnByaW50KGVwZmxfdGV4dCkKYGBgCgojIGdvZQoKYGBge3J9CmdvZV90ZXh0IDwtIHN1YnNldChkYXRhX3Bvc3RzX2NvbnRlbnQsIGlkID09ICJnb2UuY3N2IikKZ29lX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoZ29lX3RleHQpCgpnb2VfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKGdvZV9jb3JwdXMsIGMoImNhbiIsICJuZXciLCAiYW1wIiwgIndpbGwiKSkKCmdvZV93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhnb2VfY2xlYW5fdGV4dCwgMTUpCgpnb2Vfd29yZF9jYXRlZ29yeSA8LSBjKCJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiRW5nYWdlbWVudCIsICJJbWFnZSIsICJFbmdhZ2VtZW50IiwgIkVkdWNhdGlvbiIpCgpnb2Vfd29yZF9wYWlyIDwtIHBhaXJpbmcoZ29lX3dvcmRfZGljdGlvbmFyeSwgZ29lX3dvcmRfY2F0ZWdvcnkpCgpwcmludChnb2Vfd29yZF9wYWlyKQoKIyBBcHBseWluZyB0aGUgY2xhc3NpZnlfdGV4dCBmdW5jdGlvbiB0byBlYWNoIHRleHQKZ29lX3RleHQgPC0gZ29lX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IGdvZV93b3JkX3BhaXIpKQoKcHJpbnQoZ29lX3RleHQpCmBgYAoKIyBoYXJ2YXJkCgpgYGB7cn0KaGFydmFyZF90ZXh0IDwtIHN1YnNldChkYXRhX3Bvc3RzX2NvbnRlbnQsIGlkID09ICJoYXJ2YXJkLmNzdiIpCmhhcnZhcmRfY29ycHVzIDwtIGNvcnB1c19tYWtlcihoYXJ2YXJkX3RleHQpCgpoYXJ2YXJkX2NsZWFuX3RleHQgPC0gY2xlYW51cChoYXJ2YXJkX2NvcnB1cywgYygibmV3IiwgImNhbiIsICJ3aWxsIiwgInN1bW1lciIsICJ5ZWFyIiwgImZpcnN0IiwgIm1heSIsICItIiwgInJlY2VudCIsICJ5ZWFycyIsICJvbmUiLCAic2FpZCIsICJ0aW1lIiwgIm1hbnkiLCAid29ybGQiLCAiY2hhbmdlIikpCgpoYXJ2YXJkX3dvcmRfZGljdGlvbmFyeSA8LSBmcmVxX3Rlcm1zKGhhcnZhcmRfY2xlYW5fdGV4dCwgNTApCgpoYXJ2YXJkX3dvcmRfY2F0ZWdvcnkgPC0gYygiSW1hZ2UiLCAiRWR1Y2F0aW9uIiwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkVuZ2FnZW1lbnQiLCAiRWR1Y2F0aW9uIiwgIlNvY2lldHkiLCAiRWR1Y2F0aW9uIiwgIlNvY2lldHkiLCAiUmVzZWFyY2giLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgTkEsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgTkEsICJSZXNlYXJjaCIsICJTb2NpZXR5IiwgIlNvY2lldHkiLCAiUmVzZWFyY2giLCAiRWR1Y2F0aW9uIikKCmhhcnZhcmRfd29yZF9wYWlyIDwtIHBhaXJpbmcoaGFydmFyZF93b3JkX2RpY3Rpb25hcnksIGhhcnZhcmRfd29yZF9jYXRlZ29yeSkKCnByaW50KGhhcnZhcmRfd29yZF9wYWlyKQoKaGFydmFyZF90ZXh0IDwtIGhhcnZhcmRfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gaGFydmFyZF93b3JkX3BhaXIpKQoKcHJpbnQoaGFydmFyZF90ZXh0KQpgYGAKCiMgbGVpY2VzdGVyCgpgYGB7cn0KbGVpY2VzdGVyX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gImxlaWNlc3Rlci5jc3YiKQpsZWljZXN0ZXJfY29ycHVzIDwtIGNvcnB1c19tYWtlcihsZWljZXN0ZXJfdGV4dCkKCmxlaWNlc3Rlcl9jbGVhbl90ZXh0IDwtIGNsZWFudXAobGVpY2VzdGVyX2NvcnB1cywgYygi8J+RiSIsICJkYXkiLCAibmV3IiwgImNsZWFyIiwgInllYXIiLCAiY2FuIiwgIndpbGwiLCAidGltZSIsICJzcGFjZSIsICJvbmUiLCAiZmlyc3QiKSkKCmxlaWNlc3Rlcl93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhsZWljZXN0ZXJfY2xlYW5fdGV4dCwgNjApCgpsZWljZXN0ZXJfd29yZF9jYXRlZ29yeSA8LSBjKCJJbWFnZSIsIE5BLCAiSW1hZ2UiLCAiRW5nYWdlbWVudCIsICJTb2NpZXR5IiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCAiUmVzZWFyY2giLCAiRW5nYWdlbWVudCIsICJFbmdhZ2VtZW50IiwgIkltYWdlIiwgIkVuZ2FnZW1lbnQiLCAiRWR1Y2F0aW9uIiwgIkVuZ2FtZW1lbnQiLCAiRW5nYWdlbWVudCIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiRW5nYWdlbWVudCIsICJFZHVjYXRpb24iLCBOQSwgIkltYWdlIikKCmxlaWNlc3Rlcl93b3JkX3BhaXIgPC0gcGFpcmluZyhsZWljZXN0ZXJfd29yZF9kaWN0aW9uYXJ5LCBsZWljZXN0ZXJfd29yZF9jYXRlZ29yeSkKCnByaW50KGxlaWNlc3Rlcl93b3JkX3BhaXIpCgpsZWljZXN0ZXJfdGV4dCA8LSBsZWljZXN0ZXJfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gbGVpY2VzdGVyX3dvcmRfcGFpcikpCgpwcmludChsZWljZXN0ZXJfdGV4dCkKYGBgCgojIG1hbmNoZXN0ZXIKCmBgYHtyfQptYW5jaGVzdGVyX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gIm1hbmNoZXN0ZXIuY3N2IikKbWFuY2hlc3Rlcl9jb3JwdXMgPC0gY29ycHVzX21ha2VyKG1hbmNoZXN0ZXJfdGV4dCkKCm1hbmNoZXN0ZXJfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKG1hbmNoZXN0ZXJfY29ycHVzLCBjKCLwn5GHIiwgImNhbiIsICJqdXN0IiwgIndpbGwiLCAiZ2V0IiwgIndlbGwiLCAib25lIiwgImhlbHAiLCAibm93IiwgIm5ldyIsICJyZWFkIiwgImNvbmdyYXR1bGF0aW9ucyIpKQoKbWFuY2hlc3Rlcl93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhtYW5jaGVzdGVyX2NsZWFuX3RleHQsIDUwKQoKbWFuY2hlc3Rlcl93b3JkX2NhdGVnb3J5IDwtIGMoTkEsICJJbWFnZSIsICJFbmdhZ2VtZW50IiwgIkltYWdlIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCAiRW5nYWdlbWVudCIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiU29jaWV0eSIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIsICJFbmdhZ2VtZW50IiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIpCgptYW5jaGVzdGVyX3dvcmRfcGFpciA8LSBwYWlyaW5nKG1hbmNoZXN0ZXJfd29yZF9kaWN0aW9uYXJ5LCBtYW5jaGVzdGVyX3dvcmRfY2F0ZWdvcnkpCgpwcmludChtYW5jaGVzdGVyX3dvcmRfcGFpcikKCm1hbmNoZXN0ZXJfdGV4dCA8LSBtYW5jaGVzdGVyX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IG1hbmNoZXN0ZXJfd29yZF9wYWlyKSkKCnByaW50KG1hbmNoZXN0ZXJfdGV4dCkKYGBgCgojIG1pdAoKYGBge3J9Cm1pdF90ZXh0IDwtIHN1YnNldChkYXRhX3Bvc3RzX2NvbnRlbnQsIGlkID09ICJtaXQuY3N2IikKbWl0X2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIobWl0X3RleHQpCgptaXRfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKG1pdF9jb3JwdXMsIGMoIm5ldyIsICJjYW4iLCAic2F5cyIsICLigJwiLCAi4oCZIiwgIuKAlCIsICIiLCAibWF5IiwgImZpcnN0IiwgIndpbGwiLCAidXNpbmciLCAid2F5IiwgIm9uZSIsICJzY2llbmNlIikpCgptaXRfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXMobWl0X2NsZWFuX3RleHQsIDQwKQoKbWl0X3dvcmRfY2F0ZWdvcnkgPC0gYygiSW1hZ2UiLCAiUmVzZWFyY2giLCBOQSwgIkltYWdlIiwgTkEsICJFZHVjYXRpb24iLCAiUmVzZWFyY2giLCBOQSwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiUmVzZWFyY2giLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJTb2NpZXR5IikKCm1pdF93b3JkX3BhaXIgPC0gcGFpcmluZyhtaXRfd29yZF9kaWN0aW9uYXJ5LCBtaXRfd29yZF9jYXRlZ29yeSkKCnByaW50KG1pdF93b3JkX3BhaXIpCgptaXRfdGV4dCA8LSBtaXRfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gbWl0X3dvcmRfcGFpcikpCgpwcmludChtaXRfdGV4dCkKYGBgCgojIHNiCgpgYGB7cn0Kc2JfdGV4dCA8LSBzdWJzZXQoZGF0YV9wb3N0c19jb250ZW50LCBpZCA9PSAic2IuY3N2IikKc2JfY29ycHVzIDwtIGNvcnB1c19tYWtlcihzYl90ZXh0KQoKc2JfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKHNiX2NvcnB1cywgYygibmV3IiwgIuKAlCIsICJ3aWxsIiwgIndlZWsiLCAiY2FuIiwgImFtcCIsICJ2aWEiLCAibm93IiwgImZ1dHVyZSIsICJmaXJzdCIpKQoKc2Jfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXMoc2JfY2xlYW5fdGV4dCwgMzApCgpzYl93b3JkX2NhdGVnb3J5IDwtIGMoIkltYWdlIiwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIkltYWdlIiwgIkltYWdlIiwgTkEsICJSZXNlYXJjaCIsICJJbWFnZSIsICJJbWFnZSIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIkVkdWNhdGlvbiIsICJJbWFnZSIsICJTb2NpZXR5IikKCnNiX3dvcmRfcGFpciA8LSBwYWlyaW5nKHNiX3dvcmRfZGljdGlvbmFyeSwgc2Jfd29yZF9jYXRlZ29yeSkKCnByaW50KHNiX3dvcmRfcGFpcikKCnNiX3RleHQgPC0gc2JfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gc2Jfd29yZF9wYWlyKSkKCnByaW50KHNiX3RleHQpCmBgYAoKIyBzdGFuZm9yZAoKYGBge3J9CnN0YW5mb3JkX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gInN0YW5mb3JkLmNzdiIpCnN0YW5mb3JkX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoc3RhbmZvcmRfdGV4dCkKCnN0YW5mb3JkX2NsZWFuX3RleHQgPC0gY2xlYW51cChzdGFuZm9yZF9jb3JwdXMsIGMoIm5ldyIsICJ3aWxsIikpCgpzdGFuZm9yZF93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhzdGFuZm9yZF9jbGVhbl90ZXh0LCAyNSkKCnN0YW5mb3JkX3dvcmRfY2F0ZWdvcnkgPC0gYygiSW1hZ2UiLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiRWR1Y2F0aW9uIiwgTkEpCgpzdGFuZm9yZF93b3JkX3BhaXIgPC0gcGFpcmluZyhzdGFuZm9yZF93b3JkX2RpY3Rpb25hcnksIHN0YW5mb3JkX3dvcmRfY2F0ZWdvcnkpCgpwcmludChzdGFuZm9yZF93b3JkX3BhaXIpCgpzdGFuZm9yZF90ZXh0IDwtIHN0YW5mb3JkX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IHN0YW5mb3JkX3dvcmRfcGFpcikpCgpwcmludChzdGFuZm9yZF90ZXh0KQpgYGAKCiMgdHJpbml0eQoKYGBge3J9CnRyaW5pdHlfdGV4dCA8LSBzdWJzZXQoZGF0YV9wb3N0c19jb250ZW50LCBpZCA9PSAidHJpbml0eS5jc3YiKQp0cmluaXR5X2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIodHJpbml0eV90ZXh0KQoKdHJpbml0eV9jbGVhbl90ZXh0IDwtIGNsZWFudXAodHJpbml0eV9jb3JwdXMsIGMoImFtcCIsICJyZWFkIiwgIm5ldyIsICJjYW4iLCAid2lsbCIsICJ3ZWVrIiwgIndvcmsiLCAiZ3JlYXQiLCAiZGF5IiwgInZpc2l0IiwgImlyaXNoIiwgImZpcnN0IiwgImxlZCIsICJjb25ncmF0dWxhdGlvbnMiKSkKCnRyaW5pdHlfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXModHJpbml0eV9jbGVhbl90ZXh0LCA0MCkKCnRyaW5pdHlfd29yZF9jYXRlZ29yeSA8LSBjKCJJbWFnZSIsICJFZHVjYXRpb24iLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiU29jaWV0eSIsICJTb2NpZXR5IiwgIkVkdWNhdGlvbiIsICJFbmdhZ2VtZW50IiwgIkVuZ2FnZW1lbnQiLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIkVkdWNhdGlvbiIsICJJbWFnZSIsICJSZXNlYXJjaCIpCgp0cmluaXR5X3dvcmRfcGFpciA8LSBwYWlyaW5nKHRyaW5pdHlfd29yZF9kaWN0aW9uYXJ5LCB0cmluaXR5X3dvcmRfY2F0ZWdvcnkpCgpwcmludCh0cmluaXR5X3dvcmRfcGFpcikKCnRyaW5pdHlfdGV4dCA8LSB0cmluaXR5X3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IHRyaW5pdHlfd29yZF9wYWlyKSkKCnByaW50KHRyaW5pdHlfdGV4dCkKYGBgCgojIHd2CgpgYGB7cn0Kd3ZfdGV4dCA8LSBzdWJzZXQoZGF0YV9wb3N0c19jb250ZW50LCBpZCA9PSAid3YuY3N2IikKd3ZfY29ycHVzIDwtIGNvcnB1c19tYWtlcih3dl90ZXh0KQoKd3ZfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKHd2X2NvcnB1cywgYygi8J+Sm/CfkpkiLCAi8J+ZjCIsICLwn5GJIiwgImhhcHB5IiwgImRheSIsICJzZWUiLCAid2VlayIsICJncmVhdCIsICJjYW4iLCAid2lsbCIsICJ3ZWxsIiwgImtub3ciLCAibmV3IiwgIm5vdyIsICJnZXQiLCAianVzdCIpKQoKd3Zfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXMod3ZfY2xlYW5fdGV4dCwgNDApCgp3dl93b3JkX2NhdGVnb3J5IDwtIGMoTkEsICJJbWFnZSIsICJJbWFnZSIsIE5BLCAiSW1hZ2UiLCAiSW1hZ2UiLCAiRW5nYWdlbWVudCIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIsIE5BKQoKd3Zfd29yZF9wYWlyIDwtIHBhaXJpbmcod3Zfd29yZF9kaWN0aW9uYXJ5LCB3dl93b3JkX2NhdGVnb3J5KQoKcHJpbnQod3Zfd29yZF9wYWlyKQoKd3ZfdGV4dCA8LSB3dl90ZXh0ICU+JQogIG11dGF0ZShjYXRlZ29yeSA9IHNhcHBseSh0ZXh0LCBjbGFzc2lmeV90ZXh0LCB3b3JkX3BhaXIgPSB3dl93b3JkX3BhaXIpKQoKcHJpbnQod3ZfdGV4dCkKYGBgCgojIHlhbGUKCmBgYHtyfQp5YWxlX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gInlhbGUuY3N2IikKeWFsZV9jb3JwdXMgPC0gY29ycHVzX21ha2VyKHlhbGVfdGV4dCkKCnlhbGVfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKHlhbGVfY29ycHVzLCBjKCJuZXciLCAi4oCUIiwgIndpbGwiLCAiY2FuIiwgIiciLCAiZmlyc3QiLCAid29yayIsICJyZWFkIiwgImhlbHAiLCAieWVhciIpKQoKeWFsZV93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyh5YWxlX2NsZWFuX3RleHQsIDcwKQoKeWFsZV93b3JkX2NhdGVnb3J5IDwtIGMoIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCBOQSwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJJbWFnZSIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCBOQSwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIlNvY2lldHkiLCAiUmVzZWFyY2giLCAiU29jaWV0eSIsICJTb2NpZXR5IikKCnlhbGVfd29yZF9wYWlyIDwtIHBhaXJpbmcoeWFsZV93b3JkX2RpY3Rpb25hcnksIHlhbGVfd29yZF9jYXRlZ29yeSkKCnByaW50KHlhbGVfd29yZF9wYWlyKQoKeWFsZV90ZXh0IDwtIHlhbGVfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0geWFsZV93b3JkX3BhaXIpKQoKcHJpbnQoeWFsZV90ZXh0KQpgYGAKCiMgU2VudGltZW50IEFuYWx5c2lzCgojIEZ1bmN0aW9uIGZvciBlbW90aW9ucwoKYGBge3J9CmVtb3Rpb25zX21ha2VyIDwtIGZ1bmN0aW9uKHRleHRzLCBoZWlfbmFtZSl7CiAgdGV4dF9lbW90aW9uIDwtIGdldF9ucmNfc2VudGltZW50KHRleHRzKQogIAogICMgUHJvcG9ydGlvbnMgZm9yIHRleHQgdmFsdWVzIAogIGVtb3Rpb25fcHJvcG9ydGlvbnMgPC0gY29sU3Vtcyhwcm9wLnRhYmxlKHRleHRfZW1vdGlvblssIDE6OF0pKQogIAogICMgVmlzdWFsaXphdGlvbiBvZiBwZXJjZW50YWdlcyBvbiBlYWNoIGVtb3Rpb24gZm91bmQgb24gcG9zdHMKICBiYXJwbG90IDwtIGJhcnBsb3QoCiAgICBzb3J0KGNvbFN1bXMocHJvcC50YWJsZSh0ZXh0X2Vtb3Rpb25bLCAxOjhdKSkpLCAKICAgIGhvcml6ID0gVFJVRSwgCiAgICBjZXgubmFtZXMgPSAwLjcsIAogICAgbGFzID0gMSwgCiAgICBtYWluID0gc3ByaW50ZigiRW1vdGlvbnMgZm91bmQgaW4gJXMncyB0ZXh0cyIsIGhlaV9uYW1lKSwgCiAgICB4bGFiPSJQZXJjZW50YWdlIiwKICAgIHhsaW0gPSBjKDAsIG1heChlbW90aW9uX3Byb3BvcnRpb25zKSAqIDEuMikKICApCiAgCiAgdGV4dCgKICAgIHggPSBzb3J0KGVtb3Rpb25fcHJvcG9ydGlvbnMpLAogICAgeSA9IGJhcnBsb3QsICAgICAgICAgICAgICAgICAgICAKICAgIGxhYmVscyA9IHNwcmludGYoIiUuMmYlJSIsIDEwMCAqIHNvcnQoZW1vdGlvbl9wcm9wb3J0aW9ucykpLAogICAgcG9zID0gNCwgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICBjZXggPSAwLjcgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICApCn0KYGBgCgojIER1a2UKCmBgYHtyfQpkdWtlX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQoZHVrZV90ZXh0JHRleHQpLCAyKQoKcHJpbnQoZHVrZV90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihkdWtlX3RleHQkdGV4dCwgImR1a2UiKQpgYGAKCiMgZXBmbAoKYGBge3J9CmVwZmxfdGV4dCRzZW50aW1lbnQgPC0gcm91bmQoZ2V0X3NlbnRpbWVudChlcGZsX3RleHQkdGV4dCksIDIpCgpwcmludChlcGZsX3RleHQpCmBgYAoKYGBge3J9CmVtb3Rpb25zX21ha2VyKGVwZmxfdGV4dCR0ZXh0LCAiZXBmbCIpCmBgYAoKIyBnb2UKCmBgYHtyfQpnb2VfdGV4dCRzZW50aW1lbnQgPC0gcm91bmQoZ2V0X3NlbnRpbWVudChnb2VfdGV4dCR0ZXh0KSwgMikKCnByaW50KGdvZV90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihnb2VfdGV4dCR0ZXh0LCAiZ29lIikKYGBgCgojIGhhcnZhcmQKCmBgYHtyfQpoYXJ2YXJkX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQoaGFydmFyZF90ZXh0JHRleHQpLCAyKQoKcHJpbnQoaGFydmFyZF90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihoYXJ2YXJkX3RleHQkdGV4dCwgImhhcnZhcmQiKQpgYGAKCiMgbGVpY2VzdGVyCgpgYGB7cn0KbGVpY2VzdGVyX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQobGVpY2VzdGVyX3RleHQkdGV4dCksIDIpCgpwcmludChsZWljZXN0ZXJfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIobGVpY2VzdGVyX3RleHQkdGV4dCwgImxlaWNlc3RlciIpCmBgYAoKIyBtYW5jaGVzdGVyCgpgYGB7cn0KbWFuY2hlc3Rlcl90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KG1hbmNoZXN0ZXJfdGV4dCR0ZXh0KSwgMikKCnByaW50KG1hbmNoZXN0ZXJfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIobWFuY2hlc3Rlcl90ZXh0JHRleHQsICJtYW5jaGVzdGVyIikKYGBgCgojIG1pdAoKYGBge3J9Cm1pdF90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KG1pdF90ZXh0JHRleHQpLCAyKQoKcHJpbnQobWl0X3RleHQpCmBgYAoKYGBge3J9CmVtb3Rpb25zX21ha2VyKG1pdF90ZXh0JHRleHQsICJtaXQiKQpgYGAKCiMgc2IKCmBgYHtyfQpzYl90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KHNiX3RleHQkdGV4dCksIDIpCgpwcmludChzYl90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihzYl90ZXh0JHRleHQsICJzYiIpCmBgYAoKIyBzdGFuZm9yZAoKYGBge3J9CnN0YW5mb3JkX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQoc3RhbmZvcmRfdGV4dCR0ZXh0KSwgMikKCnByaW50KHN0YW5mb3JkX3RleHQpCmBgYAoKYGBge3J9CmVtb3Rpb25zX21ha2VyKHN0YW5mb3JkX3RleHQkdGV4dCwgInN0YW5mb3JkIikKYGBgCgojIHRyaW5pdHkKCmBgYHtyfQp0cmluaXR5X3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQodHJpbml0eV90ZXh0JHRleHQpLCAyKQoKcHJpbnQodHJpbml0eV90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcih0cmluaXR5X3RleHQkdGV4dCwgInRyaW5pdHkiKQpgYGAKCiMgd3YKCmBgYHtyfQp3dl90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KHd2X3RleHQkdGV4dCksIDIpCgpwcmludCh3dl90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcih3dl90ZXh0JHRleHQsICJ3diIpCmBgYAoKIyB5YWxlCgpgYGB7cn0KeWFsZV90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KHlhbGVfdGV4dCR0ZXh0KSwgMikKCnByaW50KHlhbGVfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIoeWFsZV90ZXh0JHRleHQsICJ5YWxlIikKYGBgCgojIENyZWF0aW9uIG9mIGF2ZXJhZ2Ugc2VudGltZW50IGZvciBlYWNoIEhFSQoKYGBge3J9CiMgQ29tYmluZSBhbGwgdGhlIHNlbnRpbWVudCBzY29yZXMgaW50byBvbmUgZGF0YSBmcmFtZQphbGxfc2VudGltZW50cyA8LSByYmluZCgKICAgIGR1a2VfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIGVwZmxfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIGdvZV90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgaGFydmFyZF90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgbGVpY2VzdGVyX3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICBtYW5jaGVzdGVyX3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICBtaXRfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIHNiX3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICBzdGFuZm9yZF90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgdHJpbml0eV90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgd3ZfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIHlhbGVfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXQopCgojIENhbGN1bGF0ZSB0aGUgYXZlcmFnZSBzZW50aW1lbnQgZm9yIGVhY2ggSEVJCmhlaV9hdmVyYWdlX3NlbnRpbWVudHMgPC0gYWdncmVnYXRlKHNlbnRpbWVudCB+IGlkLCBkYXRhID0gYWxsX3NlbnRpbWVudHMsIEZVTiA9IG1lYW4pICU+JQogICAgcmVuYW1lKGF2ZXJhZ2Vfc2VudGltZW50ID0gc2VudGltZW50KQoKIyBQcmludCB0aGUgdGFibGUKcHJpbnQoaGVpX2F2ZXJhZ2Vfc2VudGltZW50cykKYGBgCgojIENsdXN0ZXJzCgojIEZ1bmN0aW9uIGZvciBjb3NpbmUgbWF0cml4CgpgYGB7cn0KY29zaW5lX21hdHJpeF9tYWtlciA8LSBmdW5jdGlvbih0YWJsZSl7CiAgbnVtZXJpY2FsX2NvbHVtbnMgPC0gdGFibGUgJT4lCiAgICBzZWxlY3QoLWlkKQogIAogICMgTm9ybWFsaXppbmcgY29sdW1ucwogIG5vcm1hbGl6ZWRfY29sdW1ucyA8LSBhcy5kYXRhLmZyYW1lKHNjYWxlKG51bWVyaWNhbF9jb2x1bW5zKSkKICAKICAjIFRyYW5zcG9zaW5nIHRoZSBkYXRhIHRvIGNvbXB1dGUgc2ltaWxhcml0eSBiZXR3ZWVuIGNvbHVtbnMKICB0cmFuc3Bvc2VkX2RhdGEgPC0gdChub3JtYWxpemVkX2NvbHVtbnMpCiAgCiAgIyBDb21wdXRpbmcgdGhlIGNvc2luZSBzaW1pbGFyaXR5IG1hdHJpeAogIHNpbWlsYXJpdHlfbWF0cml4IDwtIGFzLm1hdHJpeChwcm94eTo6ZGlzdCh0cmFuc3Bvc2VkX2RhdGEsIG1ldGhvZCA9ICJjb3NpbmUiKSkKICAKICAjIENvbnZlcnRpbmcgZGlzdGFuY2UgdG8gc2ltaWxhcml0eQogIHNpbWlsYXJpdHlfbWF0cml4IDwtIDEgLSBzaW1pbGFyaXR5X21hdHJpeAogIAogICMgUHJpbnQgdGhlIHNpbWlsYXJpdHkgbWF0cml4CiAgcHJpbnQoc2ltaWxhcml0eV9tYXRyaXgpCn0gIApgYGAKCiMgQ3JlYXRpbmcgdGFibGUgZm9yIGNsdXN0ZXIgYWxnb3JpdGhtcwoKYGBge3J9CiMgSm9pbmluZyBhdHRyaWJ1dGUgcGVyY2VudGFnZV90d2VldHMgKHBlcmNlbnRhZ2Ugb2YgdHdlZXRzIG91dCBvZiBhbGwgcG9zdHMpIGFuZCBwZXJjZW50YWdlX3JlcGxpZXMgKHBlcmNlbnRhZ2Ugb2YgcmVwbGllcyBvdXQgb2YgYWxsIHBvc3RzKSBmcm9tIG51bWJlcl9wb3N0cyBhbHNvIGFkZGluZyB1bmlxdWVfaGFzaHRhZ3MgKG51bWJlciBvZiB1bmlxdWUgaGFzaHRhZ3MpIGFuZCBoYXNodGFnX3BlcmNlbnRhZ2UgKHBlcmNlbnRhZ2Ugb2YgcG9zdHMgdGhhdCBjb250YWluIGEgaGFzaHRhZykgZnJvbSBoYXNodGFncywgcGVyIEhFSQpjbHVzdGVyX3RhYmxlIDwtIG1lcmdlKHNlbGVjdChoYXNodGFncywgaWQsIHVuaXF1ZV9oYXNodGFncywgaGFzaHRhZ19wZXJjZW50YWdlKSwgc2VsZWN0KGRhdGFfcmF0aW8sIGlkLCBwZXJjZW50YWdlX3R3ZWV0cywgcGVyY2VudGFnZV9yZXBsaWVzKSwgYnkgPSAiaWQiLCBhbGw9VFJVRSkKCiMgSm9pbmluZyBhdHRyaWJ1dGUgYXZnX3Bvc3RzX3Blcl9kYXlzIChhdmVyYWdlIG9mIHBvc3RzIHBlciBkYXkpIGZyb20gcG9zdHNfcGVyX2RheSBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KHBvc3RzX3Blcl9kYXksIGlkLCBhdmdfcG9zdHNfcGVyX2RheXMpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSBhdmdfcG9zdHNfcGVyX3dlZWtzIChhdmVyYWdlIG9mIHBvc3RzIHBlciB3ZWVrKSBmcm9tIHBvc3RzX3Blcl93ZWVrIHBlciBIRUkKY2x1c3Rlcl90YWJsZSA8LSBtZXJnZShjbHVzdGVyX3RhYmxlLCBzZWxlY3QocG9zdHNfcGVyX3dlZWssIGlkLCBhdmdfcG9zdHNfcGVyX3dlZWtzKSwgYnkgPSAiaWQiLCBhbGw9VFJVRSkKCiMgSm9pbmluZyBhdHRyaWJ1dGUgYXZnX3Bvc3RzX2luX2FjYWRlbWljX3RpbWUgKGF2ZXJhZ2Ugb2YgcG9zdHMgZHVyaW5nIGFjYWRlbWljIHRpbWUpICBmcm9tIGRhdGFfcG9zdHNfYWNhZGVtaWMgcGVyIEhFSQpjbHVzdGVyX3RhYmxlIDwtIG1lcmdlKGNsdXN0ZXJfdGFibGUsIHNlbGVjdChkYXRhX3Bvc3RzX2FjYWRlbWljLCBpZCwgYXZnX3Bvc3RzX2luX2FjYWRlbWljX3RpbWUpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSBhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSAoYXZlcmFnZSBvZiBwb3N0cyBkdXJpbmcgdmFjYXRpb24gdGltZSkgZnJvbSBkYXRhX3Bvc3RzX3ZhY2F0aW9ucyBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KGRhdGFfcG9zdHNfdmFjYXRpb25zLCBpZCwgYXZnX3Bvc3RzX2luX3ZhY2F0aW9uX3RpbWUpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSB0aW1lX29mX2RheV92YWx1ZSAobnVtZXJpY2FsIHZhbHVlIHJlZmVycmluZyB0byB0aW1lIG9mIGRheSB3aGVyZSBldmVyeSBIRUkgbWFkZSBtb3JlIHBvc3RzKSBmcm9tIGZhdm91cml0ZV9ob3VyX2hlaSBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KGZhdm91cml0ZV9ob3VyX2hlaSwgaWQsIHRpbWVfb2ZfZGF5X3ZhbHVlKSwgYnkgPSAiaWQiLCBhbGw9VFJVRSkKCiMgSm9pbmluZyBhdHRyaWJ1dGUgdXJsX3BlcmNlbnRhZ2UgKHBlcmNlbnRhZ2Ugb2YgcG9zdHMgdGhhdCBjb250YWluIGFuIHVybCkgZnJvbSB1cmxfdXNhZ2UgcGVyIEhFSQpjbHVzdGVyX3RhYmxlIDwtIG1lcmdlKGNsdXN0ZXJfdGFibGUsIHNlbGVjdCh1cmxfdXNhZ2UsIGlkLCB1cmxfcGVyY2VudGFnZSksIGJ5ID0gImlkIiwgYWxsPVRSVUUpCgojIEpvaW5pbmcgYXR0cmlidXRlIGF2ZXJhZ2VfbnVtX3dvcmRzIChhdmVyYWdlIG51bWJlciBvZiB3b3JkcyBpbiB0aGUgcG9zdHMpIGZyb20gZGF0YV9wb3N0c19jb250ZW50X21ldHJpY3MgcGVyIEhFSQpjbHVzdGVyX3RhYmxlIDwtIG1lcmdlKGNsdXN0ZXJfdGFibGUsIHNlbGVjdChkYXRhX3Bvc3RzX2NvbnRlbnRfbWV0cmljcywgaWQsIGF2ZXJhZ2VfbnVtX3dvcmRzKSwgYnkgPSAiaWQiLCBhbGw9VFJVRSkKCiMgSm9pbmluZyBhdHRyaWJ1dGUgYXZlcmFnZV9zZW50aW1lbnQgKGF2ZXJhZ2Ugc2VudGltZW50IG9mIHBvc3RzKSBmcm9tIGhlaV9hdmVyYWdlX3NlbnRpbWVudHMgcGVyIEhFSQpjbHVzdGVyX3RhYmxlIDwtIG1lcmdlKGNsdXN0ZXJfdGFibGUsIHNlbGVjdChoZWlfYXZlcmFnZV9zZW50aW1lbnRzLCBpZCwgYXZlcmFnZV9zZW50aW1lbnQpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKcHJpbnQoY2x1c3Rlcl90YWJsZSkKYGBgCgpgYGB7cn0KY29zaW5lX21hdHJpeF9tYWtlcihjbHVzdGVyX3RhYmxlKQpgYGAKCiMgQmFzZWQgb24gdGhpcyBjb3NpbmUgbWF0cml4IHdlIGRlY2lkZWQgdG8gcmVtb3ZlIGF2ZXJhZ2VfbnVtX3dvcmRzIChoaWdoIHNpbWlsYXJpdHkgdG8gdW5pcXVlX2hhc2h0YWdzLCBoYXNodGFnX3BlcmNlbnRhZ2UgYW5kIHVybF9wZXJjZW50YWdlKSwgdW5pcXVlX2hhc2h0YWdzIChoaWdoIHNpbWlsYXJpdHkgdG8gaGFzaHRhZ19wZXJjZW50YWdlIGFuZCB0aGUgZGlzdGFuY2VzIG9mIGhhc2h0YWdfcGVyY2VudGFnZSB0byB0aGUgb3RoZXIgY29sdW1ucyBpcyBsb3dlcikgYW5kIGF2Z19wb3N0c19wZXJfd2Vla3MsIGF2Z19wb3N0c19pbl9hY2FkZW1pY190aW1lLCBhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSBhcmUgYWxzbyBhbGwgcmVtb3ZlZCAoZHVlIHRvIHRoZWlyIHNpbWlsYXJpdHkgdG8gYXZnX3Bvc3RzX3Blcl9kYXlzKQoKYGBge3J9CmNsdXN0ZXJfdGFibGUgPC0gY2x1c3Rlcl90YWJsZSAlPiUKICBzZWxlY3QoLWF2ZXJhZ2VfbnVtX3dvcmRzLCAtdW5pcXVlX2hhc2h0YWdzLCAtYXZnX3Bvc3RzX3Blcl93ZWVrcywgLWF2Z19wb3N0c19pbl9hY2FkZW1pY190aW1lLCAtYXZnX3Bvc3RzX2luX3ZhY2F0aW9uX3RpbWUpCgpjb3NpbmVfbWF0cml4X21ha2VyKGNsdXN0ZXJfdGFibGUpCmBgYAoKIyBGdW5jdGlvbiBmb3IgY2x1c3RlciBtZXRob2QKCmBgYHtyfQpjbHVzdGVyX21ha2VyIDwtIGZ1bmN0aW9uKG51bV9jbHVzdGVycywgdGFibGUpewogICMgRXhjbHVkaW5nIGlkIGNvbHVtbiBmb3IgY2x1c3RlcmluZwogIGNsdXN0ZXJfZGF0YSA8LSBzZWxlY3QodGFibGUsIC1pZCkKICAgIAogICMgU2NhbGluZyB0aGUgZGF0YSBmb3Iga21lYW5zIG1ldGhvZAogIHNjYWxlZF9kYXRhIDwtIHNjYWxlKGNsdXN0ZXJfZGF0YSkKICAKICBrbWVhbnNfbW9kZWwgPC0ga21lYW5zKHNjYWxlZF9kYXRhLCBjZW50ZXJzID0gbnVtX2NsdXN0ZXJzLCBuc3RhcnQgPSAxMCkKCiAgIyBFeHRyYWN0IGNsdXN0ZXIgYXNzaWdubWVudHMKICBjbHVzdGVyX2Fzc2lnbm1lbnRzIDwtIGttZWFuc19tb2RlbCRjbHVzdGVyCiAgCiAgIyBDcmVhdGUgYSBkYXRhIGZyYW1lIGNvbWJpbmluZyBvcmlnaW5hbCBkYXRhIHdpdGggY2x1c3RlciBhc3NpZ25tZW50cwogIGNsdXN0ZXJlZF9kYXRhIDwtIGNiaW5kKGNsdXN0ZXJfdGFibGUkaWQsIGNsdXN0ZXJfZGF0YSwgY2x1c3RlciA9IGNsdXN0ZXJfYXNzaWdubWVudHMpCiAgCiAgY2x1c3RlcmVkX2RhdGEgPC0gY2x1c3RlcmVkX2RhdGFbLCBjKCJjbHVzdGVyX3RhYmxlJGlkIiwgImNsdXN0ZXIiKV0KICAKICBwcmludChjbHVzdGVyZWRfZGF0YSkKfQpgYGAKCiMgRnVuY3Rpb24gdG8gZGlzY292ZXIgYmVzdCBudW1iZXIgb2YgY2x1c3RlcnMKCmBgYHtyfQplbGJvd19tYWtlciA8LSBmdW5jdGlvbih0YWJsZSl7CiAgY2x1c3Rlcl9kYXRhIDwtIHNlbGVjdCh0YWJsZSwgLWlkKQogIHNjYWxlZF9kYXRhIDwtIHNjYWxlKGNsdXN0ZXJfZGF0YSkKICAKICB3c3MgPC0gdmVjdG9yKCkKICByYW5nZSA8LSAxOjEwCiAgCiAgZm9yIChrIGluIHJhbmdlKSB7CiAgICBrbWVhbnNfbW9kZWwgPC0ga21lYW5zKHNjYWxlZF9kYXRhLCBjZW50ZXJzID0gaywgbnN0YXJ0ID0gMTApCiAgICB3c3Nba10gPC0ga21lYW5zX21vZGVsJHRvdC53aXRoaW5zcwogIH0KICAKICBlbGJvd19kZiA8LSBkYXRhLmZyYW1lKGsgPSByYW5nZSwgV1NTID0gd3NzKQogIGdncGxvdChlbGJvd19kZiwgYWVzKHggPSBrLCB5ID0gV1NTKSkgKwogICAgZ2VvbV9saW5lKCkgKwogICAgZ2VvbV9wb2ludCgpICsKICAgIGxhYnMoeCA9ICJOdW1iZXIgb2YgQ2x1c3RlcnMiLCB5ID0gIldpdGhpbi1DbHVzdGVyIFN1bSBvZiBTcXVhcmVzIChXQ1NTKSIsCiAgICAgICAgIHRpdGxlID0gIkVsYm93IE1ldGhvZCBmb3IgT3B0aW1hbCBrIikgKwogICAgdGhlbWVfbWluaW1hbCgpCn0KYGBgCgojIFBsb3Qgb2YgRWxib3cgTWV0aG9kIGFuZCBzZWxlY3Rpb24gb2YgYmVzdCBudW1iZXIgb2YgY2x1c3RlciB0byB2aWV3IGhvdyBIRUlzIGFyZSBncm91cGVkCgpgYGB7cn0KZWxib3dfbWFrZXIoY2x1c3Rlcl90YWJsZSkKY2x1c3Rlcl9tYWtlcig0LCBjbHVzdGVyX3RhYmxlKQpgYGA=